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ABSTRACT (SNSI) 



1. Summary 

As a follow-on to the mission analysis study of the 
inspection of orbiting satellites suspected of carrying nu¬ 
clear weapons (SAMSO TASK 11), Los Alamos Scientific 
Laboratory (LASL) conducted a detailed experimental 
and theoretical investigation of nuclear sensor techniques. 
The objectives were to (1) define for the 1970-1980 peri-, 
od the observables constituting weapon detection and the 
effects of potential countermeasures, (2) analyze the sen¬ 
sors recommended in Task 11 and explore new sensor 
concepts, (3) select preferred techniques, and (4) define 
the performance of the selected techniques as well as the 
imposed mission constraints. 

Unique identification of weapon fission fuels and 
discrimination from other radioactive nuterials is 


accomplished primarily by high-resolution [Ge(Li)] 
gamma spectrometry. High-resolution spectra from eight 
US weapons representing the yield range from 50 kt to 
10 Mt were measured and analyzed for characteristic ener¬ 
gies and for absolute flux as a function of distance and 
view angle in order to specify the spectrometer require¬ 
ments. Low-resolution gamma spectra and flux from a 
number of other weapons have also been evaluated. Simi¬ 
lar data for neutron flux were obtained in a few cases 
because of the redundancy that can be provided by neu¬ 
tron deteaors in the identification of material and estima¬ 
tion of quantities. 

Unique identification of deuterium as the principal 
thermonuclear fuel requires photo interrogation at ener¬ 
gies (typically 2.62-McV gamma or 4-McV brems- 
strahlung) above the D(7,n) threshold, with subsequent 
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neutron detection. Discrimination from interfering 
Be(7,n), which has a lower threshold energy, is made 
possible by changing source energy (to 1.69-MeV gamma 
or 2.2-MeV bremsstrahlung). by determining neutron end¬ 
point energy with ^He neutron spectrometry, or by rn'.as- 
uring neutron average energy with moderated ^He neu¬ 
tron detectors. With these incident gamma energies, which 
arc below the photofission threshold, the same neutron 
spectral techniques allow determination of the coupling 
between fission and fusion fuels through the D(7,n) plus 
U(n,0 reactions. I ^ ^ 


_ TTic ^mma sources for the experimental 

photo interrogation studies were either selected radio¬ 
isotopes or one of several pulsed bremsstrahlung sources, 
and an extensive series of small sample, mockup, and 
weapon irradiations were done to establish neutron yields 
and spectral characteristics. Neutron spectrometer devel¬ 
opment was also necessary to establish instrumentation 
limitations. The final calculational and experimental re¬ 
sults from this sequence of investigations are presented 
for six typical US weapons and USSR weapon models in 
the 1- to 25-Mt range. 

The identification of weapons by the characteristic 
arrangement of essential components is accomplished pri¬ 
marily by radiographic imaging. Essentially the same 
sources (energy and intensity) appropriate to photointer¬ 
rogation are applicable to radiography, although the 
source and detector (radiographic screen) must now be on 
opposit e sides of the target. T 






\ AjTcv^lua- 

don o{ radio^phic detector components was conducted 
first in order to select equipment for further paranKtric 
studies on source and geometry requirements. Seven rep¬ 
resentative weapons were used as targets in the course of 
this work, which involved a wide variety of isotopic and 
pulsed sources, j j 


L 


These essential techniques are further fortified by 
the use of pulsed sources (7.5-MeV bremsstrahlung or 
14-MeV neutrons) for inducing fissions in the target weap¬ 
on. The fission neutron yield and/or time history provides 
another unique identification of fissionable materials and 
in some cases gives significant design information. The 
availability of 7.5-MeV bremsstrahlung also extends the 
radiographic capability to heavier targets, provides some 


limited material identification through characteristic 
gamma decay of isomeric states, and, in conjunction with 
neutron spe<r.rometry, may help confirm thermonuclear 
design inforr Jition. Many experiments and calculations 
were conducted in these general areas to assist in defining 
sources, instrumentation, and operational requirements. 

Because there is no optimum or preferred system 
but rather a continuum encompassing many variables in 
areas such as weapon design, space background, and oper¬ 
ational constraints, the results are presented in four dis¬ 
crete steps. The simplest sensor package (SP-1) which 
provides some capability for each of the principal tech¬ 
niques mentioned above (gamma spearometry, radio- 
graphic imaging, and photointerrogation above the deute¬ 
rium photodisintegratio n threshold) is specified a s^jthc 
minimu m useful sptemJ __ 




SP-4, the most complex sensor pacl^ge to be consISored^ 
at all practical (including the techniques above plus, pho- 
cofission yield and time-history measurements),^ 

r ' -- -■ hO) 


I 


_mo^rcaiistic package 

for actual hardware development probably lies somewhere 
between these extremes. 

This brief summary of the recommended tech¬ 
niques, and possible options in sensor hardware to imple¬ 
ment the techniques with varying cost and confidence, 
indicates satisfactory attainment of the study objectives. 
General conclusions and recommendations regarding sub¬ 
sequent efforts can be found in Sec. VII. It is important 
to emphasize that there is no optimum system, nor will 
any of the systems described here result in complete 
confidence in weapon detection under all postulated con¬ 
ditions. The possible sensor systems represent a con¬ 
tinuum, and those presented in this report indicate a best 
estimate of the practical limits to that continuum. 


II. Introduction 

During the period from June 1970 through Febru¬ 
ary 1971, Task 11 was conducted by Space and Missile 
Systems Organization (SAMSO) using two industrial con¬ 
tractors, with technical direction and ancillary supp ort 
from Aerospace Corporation, 


!! 
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Ubontories provided advice and assistance to SAMSO 
before the contraaors were seleaed and also during this 
mission analysis effort, in particular with regard to the 
nuclear sensor and interrogation subsystems and to basic 
nuclear weapon data. Predictably, one of the study find¬ 
ings emphasized the need for a follow-on nuclear sensor 
analysis and experimenudon effort to clarify and quan- 
dze system parameters in the nuclear package. Thus, in 
March 1971 a nine-month theoretical and experimental 
program was undertaken by LASL. Progress during this 
study was reported periodically in a working group com¬ 
posed of representatives from LASL, SAMSO, Aerospace, 
Division of Military Application (DMA), AEG 
Albuquerque Operations Office (ALO), Lawrence Liver¬ 
more Laboratory (LLL), and Sandia Laboratories. The 
present repon constitutes the final summary and outlines 
the conclusions of the study. 

The statement of work agreed upon at the outset 
listed four separate technical areas for detailed investiga- 
don and provided guidance on the relative effort to be 
expended on each. Briefly, these areas include: 

(1) (10%) Defiiiidon of weapon design options dur¬ 
ing the 1970-1980 period, the possibility of tai¬ 
loring weapon designs to defeat detection and 
identification by the proposed nuclear sensors, 
and the potential penalties associated with such 
countermeasures. 



with cenain infornution to identify further their major 
characteristics. These will be referred to without further 
amplirication in later sections of this report. Most of the 
small sample and weapon component work was carried 
out at these LASL facilities. However, to get additional 
passive and active data on a variety of weapons conven¬ 
iently, several field trips were made to the AEG fabri¬ 
cation and assembly plant at Amarillo, Texas (Pantex). 
The linac used routinely there for inspection and quality 
control was used as a bremsstrahlung source and is also 
listed in Table 1. 

The complex nature of the interface between the 
nuclear instrumentation used to detect and identify weap¬ 
ons and the spacecraft required to carry that instru¬ 
mentation precluded any real sensor choices at this stage 
in the overall development program. Rather, the wide 
latitude in system synthesis has been indicated by illus¬ 
trating several potential systems with increasing levels of 
complexity,* size, and information return. The possible 
sensor systems really represent a continuum, and those 
presented in this report indicate a best estimate of the 
practical limits to that continuum as well as several inter¬ 
mediate examples. The correct inference is that hard 
engineering choices would have to be made by negotiation 
between the various requirements on a spacecraft-booster 
combination, and that a flyable sensor system would 
probably lie somewhere between the minimum-useful and 
maximum-complexity cases illustrated here. 


(2) (15%) Tabulation of gamma spectra and flux and 
neutron flux for representative weapon designs, 
including distance and weapon aspect angle con¬ 
siderations. 

(3) (45%) Experimental investigation of radiographic 
techniques for this application using a large scin¬ 
tillating screen, a wide-angle lens, high-gain image 
intensifier, and photographic recording. The pa¬ 
rameter study was to emphasize the minimum 
radiation intensity and energy required for posi¬ 
tive detection. 

(4) (30%) Application of a variety of neutron and 
photo interrogation techniques with both experi¬ 
mental and theoretical delineation of source and 
detector requirements and with typical weapon 
responses in neutron and gamma yields and time 
histories. 

Because the experimental study was short-term and 
funded at a modest level ($250K), no new source facilities 
were contemplated for the aaive interrogation studies. 
Rather, a variety of existing installations at LASL were 
used without modification. Table 1 lists the sources along 


III. Weapon Designs 

The general design features and characteristics that 
may affea nuclear detection and identification systems 
for present weapon technology have been described else¬ 
where^ and will not be reproduced here. However, the 
specific strategic warheads (RV’s) which were in the cur¬ 
rent US inventory and which were used in the experi¬ 
ments conduaed during the study are typical pegpoints in 
a wide range of generalized designs (50 kt to 10 Mt). The 
relevant parameters for these are summarized in Table 11 
and will also be referred to in various later sections of the 
report. The weapons presented in Table II actually span a 
decade or more of modern developments and certainly 
represent the current state of the technology. Table II 
also lists characteristics as calculated for one plausible set ^ 
of warheads matching thc | ^pabilities. These designs ^ 
arc used for reference points in later sections where sensor 
performance is calculated. 

In addition to the stockpile weapons and extra¬ 
polated designs, some test devices and/or subassemblies 
arc referenced in specific experiments and calculations. 

The significant characteristics of these will be given in the 
appropriate sections. 
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TABLE I 

ACTIVE INTERROGATION FACILITIES 



Machine 

Beam Energy 
and Type 

Target and 
Thickness 

Spot Size 
(in.) 

Pulse Width 

rep Rate 

SCAPP 

GMX-11 flash x-ray 
Single Cavity 

PHERMEX Prototype 

3-7 MeV 
bremsstrahlung 

0.010-0.043-in. 

W 

1/8-1 

0.2 Msec 

Single pulse 
to 1 pps 

Febetron 

GMX-1 flash x-ray 

2.3 MeV 
bremsstrahlung 

0.025-in. Ta 

1/4 

0.020 Msec 

Single pulse 

Million 

Volt 

Machine 

GMX-1 

1 MeV 

bremsstrahlung 

0.002-in. W 

3/8 

Sinusoidal 
^100 Msec 

1000 cps 

Betatron 

GMX-1 

22 MeV 
bremsstrahlung 

Internal 

^.100-in. Pt 

0.008 

1 Msec 

180 cps 

Pantex 

Linac 

Linear accelerator 

7-10 MeV 
bremsstrahlung 

0.040-in. Pt 

1/16 

^ 5 Msec 

15-250 cps 

EPA 

MP-3 EPA 

4-27 MeV 
bremsstrahlung 

0.002-in. W 

<0.4 

50-500 Msec 

120 cps 

Cockcroft- 

Walton 

Deuterium 

accelerator 

350 keV 
deuterium 
(14-MeV neutrons) 

Deuterium, 

Tritium 

NA 

1 nsec-dc 

10* cps 

Zipper 

Sandia high-output 
neutron generator 

180 keV 
tritium 

Deuterium 

NA 

5-30 Msec 

1-10 cps 


(14-McV neutrons) 




As presently foreseen, advanced development over 
the next decade is directed along two principal lines 
insofar as strategic weapons arc concerned. The first, and 
most predictable, effort concerns the refinement of pri¬ 
maries to smaller sizes, thus permitting small improve¬ 
ments in yield-to-weight ratios at the low end of the 
strategically significant yield range. Because of the large 
leverage which results from this factor, particularly in 
(systems, this approach is expected to be 
continued for at least as long as nuclear testing contmues. 
The nature of these dev elopm ents, 




IS suen as to 


have relatively little effect on nuclear detection and iden¬ 
tification techniques. Modest increases in amounts of 
fissile material present would make detection easier with 


some of the proposed techniques. On the other hand, 
increased emphasis on hardening against ABM tactics 
^ould somewhat degrade other sensor methods^ 


r: 




\ 


Application to strategic warheads appears unlikely. 
In any event, secondaries are not expected to change 
radically from those presently stockpiled in the strategic 
systems,] - _ 

f — ^ . •• 

\ Thus, this line of future dcvel- 

"dpmchtTS'hot expected To produce any large surprises in 
j^e detection and identification game in the 197P’s.^-- 
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WEAPON DESIGN CHARACTERISTICS 
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The moti^don for the strategic designer perhaps lies in 
an increase by a factor of three or four in yield-to-weight 
rados, but the probability of successful deployment of 
such a weapon in the 1970-1980 period is low and unpre¬ 
dictable. Because of this uncertainty, the design of an 
inspection system should not be influenced strongly by 
characterisdes ^ 

Assuming typical and current weapon technology, 
the important criteria for detection and identificadon 
can be listed. While the order of listing is not necessarily 
significant, two categories are used which to some extent 
separate the deteaion and idendfication criteria from 
those of acquiring significant additional intelligence data. 
Obviously there is considerable overlap in these two cate¬ 
gories as well as in the individual criteria. Some items as 
listed are not entirely separable except under the limita- 
dons of specific techniques, nor are all items relevant to 
all potential warheads. However, all are discernible to 
some degree by the inspection methods proposed in this 
report. 

(1) (a) Unique identificadon of the essential fission 
and fusion fuels: ^Pu, and deute¬ 

rium. 


(0 Measurement of nuclear material isotopic 
rad os for production cycle information. 

(g) Determinadon of reaedvity of subcritical 
configurations by dynamic response. 

These criteria, if completely determinable, would 
give high confidence to any statement as to the presence 
of nuclear warheads and would with reasonable confi¬ 
dence add to the completeness and accniracy of current 
US estimates of foreign weapon technology. The com¬ 
plexity of the list requires a combination of techniques to 
maximize the chance of obtaining reliable data on any 
given target. As will be seen later in the correlation of 
sensor packages to detection criteria, the overall success 
rate is related directly to the number and sophistication 
of the detecdon techniques employed and to the relative 
absence of countermeasures. 


IV. Countermeasures 




\ 


(b) Determinadon of approximate quantity and 
distribution of fissile and fusion fuels. 


f 


(c) 


I 

(2) (a) Idendfication and configuration of high-Z^ 
materials in nuclear components and radiadon* 
case. 

(b) Determinadon of radiadon channel, HE, and 
RV configuradons. 

(c) Measurement of approximate attenuation o£ 
system in unique locadons. 

(d) Identification of RV- and warhead*i 
hardening techniques. 

( 

(e) Determination of significant fuzing and fir-j 
ing details that might be associated with boosdng 
or safing components. 
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TABLE 111 

DESIGN-INDUCED COUNTERMEASURES 


Approximate Yield Penalty 
for Warhead of Size: 


Affeaed Inspection Countermeasure 

Technique _ Effectiveness 100-500 kt 0.5-5 Mt > 5 Ml 


hU) 


I 



I 


h^) 

_ _ _i 

A useful reference exhibiting spectra of many of the 
above-mentioned isotopes is Idaho Nuclear Corp. report 
lN-1448 Rev. (January 1971). 


2. Gamma Spectra of US Weapons, The following 
figures display the basic features of the gamma spectra 
from several types of US weapons. Some of these data 
have been published previously,’ but separate spectra 
taken at the primary and secondary arc shown here to 
illustrate the relative intensity of the various lines. These 
particular spectra arc all taken close to the weapon 
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TABLE IV 



ADD-ON COUNTERMEASURES 


Approximate Weight Penalty 
for Warhead of Size: 


Countermeasure 

Type 


Affected Inspecdon 
Technique 


Countermeasure 

Effeedveness 100-500 kt 0.5-5 Mt > 5 Mt 
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3 ft) and cannot be used to infer count rates at larger 
distances. However, an excellent idea of the spectral fea¬ 
tures observed at large distances can be obtained by 
combining visually the primary and secondary spectra. 
Table II gives the basic design features of these devices. 

Figure 1 shows the spectrum from the W-62 with 
the detector looking in the channel between the primary 
and secondary, in this and subsequent spectra, unless 
otherwise stated, the 6% Ge(Li) detector (relative to 3- by 
3-in. Nal at 1.33 MeV and 25 cm distance) is shielded and 
collimated to a 74® field of view (FOV). This spectrum of _ 
the W-62 is the same as that seen at a distance) 








Spec tra from th e Mk -28 are exhibited in Fig. 2.^ 




The Mk-43 spectra in F^. 3 are similar to those of 
the Mk-28.1 ' , 




1 \ 
r— 






Similar spectra arc observed from the 

W-47 and W-58. 

Spectra from the pr imary and second ary of the_ 
Ai k-53 are sh own in Fig. 5.^ 

I I 


(L)i ^ 




In Fig. 6 spectra from the primary (top), and sec¬ 
ondary (bottom) of the W-59 arc shown. These spectra 
were taken with a 41® FOV, shielded, 4.5% relative effi¬ 
ciency Ge(Li) detector 20 in. from the axis of the weap¬ 
on. Counting times wcr^4000 sec, and the weapon was in 
1^ its shiODins container.. ( 

i 

These figures emphasize the design-dependent range 
of gamma outputs. As is well known, there is no one basic 
design parameter (size, weight, yield, etc.) that can be 
correlated directly with the gamma output. For examplCj^^ 
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FSfr4. the spectrum from the priimry of the 
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Fig. 1. 


W-62 gamma spectrum, 32 in. from axis between j 
primary and secondary, 1000 sec. 
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Fig. 2. 

Mk-28gamma spectra, 28 in. from axis, 1000 sec. 
Top plot is at primary, bottom at secondary. 
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Mk-43gamma spectra, 27 in. from axis, WOO sec. 
Top plot is at primary, bottom at secondary. 
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Ftg. 4. 

Vi-SSgamma spectra, 44 in. from axis, 1000 sec. 
Top plot is at primary, bottom at secondary. 
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Fig. 5. 

Mk-53 gamma spectra, 30 in. from axis, WOO sec. 
Top plot is at primary, bottom at secondary. 
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W-59gamma spectra, 20 in. from axis, 4000 sec with 4.5% Ge(Li). 
Top plot is at primary, bottom at secondary. 

















^ata to be shown later will also illustrate the strong 


view-angle dependence of these spectra. 

These spectra indicate that the lines that can be 
used for isotopic identification in practical cases are the 
characteristic plutonium and lines. The 2614-keV 
line is there, but its presence, in general, cannot be unam¬ 
biguously assigned to a specific isotope, because it can be 



observed from 93.5% “’Pu, and ^^Th. 

rr-^-^- 



Fig. 7. 


W-59 pbotopeak flux vs distance, plutonium photo- 
peaks, Nal and Ge(Li) detectors. 


3. Experimental Results. The major experimental 


gamma-spectroscopy effort in this program was devoted 
to measuring photopeak fluxes from various US weapons 

systems as a function of distance and view angle. A between the side and the primary end of the device. A 


discussion of some characteristics of Ge(Li) detectors is similar comparison for the D-38 lines for side incidence is 
relevant to the effort and is given separately in Appendix shown in Fig. 8. 

A. The measurements were made with 6% Ge(Li) and 3- The scatter of the points is indicative of the magni- 



by 3-in. Nal detectors, and net counts in each photo peak tude of errors that arise in the measurement. 


were obtained by subtracting a linear extrapolation of the 
Compton continuum from underneath the photopeak. 
The efficiency curves of the appendix were then used to 
convert net photopieak counts to photopeak flux. 


For the Nal data, the 400-keV ^Pu complex 
appears as one peak. The flux was computed for the 
entire complex using the detector efficiency value at 
400 keV. The high resolution data of Qine^ shows that j 
the 375- and 414-keV thin sample intensities are equal. 

The sum of the 375- and 414-keV intensities in turn is 
equal to approximately one-half the total intensity be- | 
tween 300 and 450 keV, the range covered in the Nal i 
peak. With these arguments for a basis, the 375- and i 
414-keV photopeak fluxes or counting rates were approx¬ 
imated by taking each to be equal to one-founh of the ‘--- 

total flux in the 400-keV complex. This approximation 
neglects minor differences in anenuation in the weapon; 
however, these differences tend to be cancelled in the 
detector counting rate because of the slope of the detec¬ 
tor efficiency curve. 

Data were taken at LASL on the W-59 in its ship¬ 
ping container with both Nal and Ge(Li) detectors in a 
33^ FOV geometry. Measurements as a function of dis¬ 
tance out to about 30 ft were nude at four different view 

angles. For two of the view angles both Ge(Li) and Nal |i 

data were taken. In Fig. 7, a comparison of the Ge(Li) \ _ 

and Nal data is shown for the 375- and 414-kcV lines as Fig. 8. 

well as for the entire 400-keV complex. The view angle W-59 pbotopeak flux vs distance, pbotopeaks, 

here is looking at the center of the channel at a 45*’ angle Nal and GeiLi) detectors. 
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TABLE V 


PHOTOPEAK FLUX VS VIEW ANGLE, W-59 AT 20 FT 



Flux ( 7 /cin^/scc) 


Fig. 10. 

Mk’43 pbotopeak flux vs distance, and 

photopeaks, Nal detector. 




These studies were carried out on four different 
designs (W-62. W- 68 , Mk-5 3, Mk-43) at Pantex and Figs. 9 
through 11 show the results of some of the measure* 
ments, all taken with a 3- by 3-in. Nal detector in 51® 
FOV geometry. Typical spectra from which these data 
TOin t s were extracted are shown in Figs. 12 thr ough 16. 


These illustrated spectra were all taken at 20 ft for count 
times of 10 to 20 min and represent realistic spectra that 
would be obtained in a low background inspection for 
reasonable counting times if a 3- by 3-in. Nal detector 
were used. The results of the flux measurements for a 
distance of 20 ft are tabulated in Table VI. 


IT 



Fig. 11. 

Mk-53 photopeak flux vs distance, photopeak, 
Sal detector, end-on view of secondary. 


I 




( 2 ^ 
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Fig. 9. 

W-62 and W-68 photopeak flux vs distance, ^^Pu 
photopeaks, Sal detector. 


Fig. 12. 

W-68 gamma spectrum, 3- hy 3'in. Sal at 20ft, 
10-min count. 
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^ Fig. 13. 

W'62 gamma spectrum. J- by 3’in. Nal at 20 ft, 
20-min count. 


Fig. 14. 

Mk’43 gamma spectrum, 3- by 34n. Nal at 20ft, 
lO-min count. 
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Fig. 15. 

Mk-53 gamma spectrum, 3- by 3^n. Nal at 20 ft, 
20-min count, end-on view of secondary. 
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Fig. 16. 

Background gamma spectrum, Nal detector, 20-min 
count. 

All systems were viewed from the side except the 
Mk-53, which was viewed from the secondary end. All 
collinudon was to 51^ FOV except for a 33^ FOV used 
for the W-59. All systems were bare except the W-59 and 
Mk-53, which were in their shipping containers. 

Questions regarding collimator FOV can be an s wer- 
ed rather simplyj[ 

r'] ! 

I ^l^liepi^hden^’will not begin 

to set in until the PbV oecomes larger than the object 
under view, and the distance from the device is several 
times the largest linear dimension of the source. Ideally, 
the collimation should be tight enough for the inspection 
distance chosen so that the object viewed nearly fills the 

TABLE VI 

MEASURED PHOTON FLUX AT 20 FT 
(7/cm ^/sec) 


Gamma Energy 
(keV) 


i(0 


Device 


400 complex 


766 


1001 


% 


^‘^Parentheses indicate a very approximate value. See dashed line 
in Fig. 10 and spectrum in Fig. 14. All further calculations using 
this number are also approximate. 
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Wum colGn^on dq>cnd^ of courscrbn the size of the 
target and the inspection distance and implies a significant 
weight associated with the shield. 

The Task 11 mission analysis included calculations 
of photopeak flux from three specific weapon models. 
None of the present measurements were made on the 
models calculated. However, from design similarity argu¬ 
ments the measured results can be roughly scaled to 
compare to the contractor calculations. Agreement by a 
factor of two to three is considered reasonable consider¬ 
ing the calculational approximations. 

4, Minimum Counting Time for Detection. From 
the measured fluxes at 20 ft (Table VI) and the 
18% Ge(Li) efficiency curve in the appendix, the photo¬ 
peak count rates expected for the various sources were 
ubulMd and are shown in Table VII. 


Two factors influence the conversion of these count 
rates to detection time. The simplest is that there are no 
universally accepted criteria for minimum detection. The 
work of Currie*® points out that differences of over a 
factor of 30 can appear for reasonable but differing cri¬ 
teria. The most serious problem, though, is the strong 
dependence of the minimum counting time on back¬ 
ground. The details of the background expected in a 
{detector are difficult to calculate because 

TABLE VII 

PHOTOPEAK COUNT RATES 
(counts/min) 

18% Ge<U)-at 20 ft 

Gamma Eneri^ 

(keV) 


of the number of parameters involved and our incomplete 
knowledge of them. The space radiation background is 
discussed further in Sec. VI and in Appendix C, but for 
the present minimum detection time calculations are 
based on the backgrounds observed in shielded detectors 
in aaual measurements on weapons. Here background is 
defined as the continuum on which the photopeaks ride. 
In the absence of large external sources, most of the 
background under the photo peaks arises from Compton 
scattering in the source and the detector. In this low 
background case both the magnitude and shape of the 
background arc source-dependent. Reduction of some of 
this background by Compton suppression techniques is 
discussed in this seaion. 

Backgrounds observed in measurements on the 
W-59 at 20 ft with a 6% Ge(Li) detector have been scaled 
up to the larger 18%Ge<Li) detector (shielded) and are 
listed in Table VIII for 2048 channels at 1.4 keV/channel. 

Sophisticated computer techniques that search spec¬ 
tra for photopeaks and then analyze the peaks for energy 
and area are commonly available. One such code** 
(SAMPO) was used in this program to develop criteria for 
minimum area photopeak detection.*^ This work consid¬ 
ered: (l)the probability that background fluctuations 
would be identified incorrectly as true photo peaks and 
(2) the probability that true photopeaks would be dis¬ 
carded in the process of discriminating against fraudulent 
peaks from background fluctuations. These probabilities 
were analyzed in terms of the background under the 
photopeak, and minimum counting time criteria were 
developed. However, this work considered only one width 
of photo peak, and the results are not universally appli¬ 
cable. The analysis is being extended, and eventually these 
computer search techniques will probably be more sensi¬ 
tive than the method outlined below. The latter, and 
more conservative, approach is adequate and perhaps de¬ 
sirable for system planning purposes. 

The detection criterion assumed here is that the 
photo peak height is five times the standard deviation of 
the background points under the peak. The peak area of 
the assumed Gaussian peaks is simply related to the full 

TABLE VIII 


Device 


375 


414 


766 


1001 


ESTIMATED MINIMUM BACKGROUND COUNT RATES 


Ey 

(keV) 


(«), 


Parentheses indicate very approximate value. 


Counts/min/channel 


375 

3 


414 755 


2.5 0.85 


1001 


0.25 
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width at half maximum (FWHM) and height of the peak 
by the equation: 

Area * 1.064 x Height x FWHM. 

Using this relation and the fact that the standard devia¬ 
tion of a background point under the peak is 
the following expression can be derived for the minimum 
counting time: 

28.30 FWHM* 

- , ( 1 ) 

background count rate under 
peak (counts/channel-min), 

photopeak full width at half 
maximum (channels), 

integrated photopeak count 
rate (counts/min). 




where R^^i^ 


FWHM 


- 


The criteria of Eq. (1) applied to the count rates of 
Table VII and the background of Table VIII produce the 
minimum detection times listed in Table IX. A FWHM of 

TABLE IX 

MINIMUM DETECTION TIME IN MINUTES 


Energy 

(kcV) 


Device 375 414 766 1001 

•T -- - -- - --- 

1 


four channels has been assumed. The larger of the follow¬ 
ing times is listed in the uble: (1) the time from Eq. (1), 
or (2) the time to collect 50 counts in the peak. The 
choice of 50 counts is somewhat arbitrary but is governed 
by the desire for reasonable statistics in determining ap¬ 
proximate peak areas. In cases where the latter is larger, 
the Eq. (1) time is listed in parentheses for possible use in 
scaling with background. Note that the minimum detec¬ 
tion time defined by Eq. (1) scales direaly with back¬ 
ground. 

These detection times are not overly conservative 
when it is considered that enough counts must be accu¬ 
mulated so that meaningful relative intensity measure¬ 
ments can be made to avoid confusion from possible 
countermeasures. Computer search techniques can pro¬ 
duce minimum count times that are lower than those in 
Table IX, but the resulting peak area accuracy docs not 
allow relative intensity measurements. It has been 
shownthat count times for quantitative intensity deter¬ 
minations can be greater than the minimum times for 
detection by as much as a factor of 10.___ 



t 


/ 


5. Compton Suppression Spectrometers. A non- 
exhausdve search of the literature on Compton suppres¬ 
sion and anticoincidence shield techniques was unsuccess¬ 
ful in the sense that very linlc information was found on 
systems that would be suitable for satellite applications. 
The majority of the literature describes massive fixed- 
installation laboratory systems. |_ _ 
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spccmim from a 0.003-in. 
chat from a] 


.[As an exampc,‘ tic 
®U foil was compared with 


% 




I Ultrahigh resolution Si(Li) and Ge(Li) x-ray detec- 

I tors with responses down to a few kiloelectron volts are 
I common laboratory items. The small size of these low- 
energy detectors precludes their uses 

at the present time, but such a detector has 
beenHsecTto illustrate the above points.^ 


U 




A mantled detector is considered in greater detail in 
Sec. VI and in Appendix C. 


6. Characteristic X-Rays. Detectable x-ray emis¬ 
sions from nuclear weapons have been studied to deter¬ 
mine the quantity and quality of intelligence information 
that might be obtained from such spectra. 

There are two principal methods of producing 
X rays in the radioactive decay of weapons materials: 


(a) Internal conversion in the daughter nucleus after 
the parent decay produces x rays characteristic of 
the daughter, i.e., uranium K x rays dominate the 
K x-ray spectrum of thin plutonium sources. 


(b) 


Fluorescent K x rays are produced in thicker 
sources from the photoelectric absorption of 




higher-energy gamnu radiation emitted by the 
same material. This effect produces K x rays 
characteristic of the material itself. For example, 
in a thick uranium source, the fluorescent urani¬ 
um X rays dominate over the thorium x rays 
expected from internal conversion in the daugh¬ 
ter nucleus. 
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Fig. 17. 

W-68 characteristic x-ray spectrum, primary end. 


B. Passive Neutron Detection. The significant neu¬ 
tron emitters among common weapons materials are “U 
and and these isotopes emit approximately 

0.014 n/g-scc and 900 n/g-sec respectively from spontane¬ 
ous fission.! 


; 

) 


*"Polycthylene-modcratcd ^Hc neutron detectors 
have been developed with intrinsic efficiencies for 1-MeV 
neutrons of nearly 20%.®**^*^ For example, a detector of 
this type (19 ®Hc tubes in a 20- by 24- by 3-in. 
polyethylene matrix) was used to detect neutrons from 


the W-59. The count rate as a function of distance for a 
side-on position is shown in Fig. 18. For this detector, the 
local background ranges from 5 to 10 cps. 

While passive neutron detection is desirable, it is 
also redundant since passive gamma techniques give a 
more positive indication of the presence of ^Pu. How¬ 
ever, passive neutron detection is provided as an option 
for little cost by using the neutron detector proposed for 


I 



W-59 neutron count rate vs distance, 20- by 24- by 
3-in. moderated ^ He detector. 
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^ Examples of the large colleaion of radiographs, 

including a variety of weapon types and a variety of 
exposure conditions, arc given in the following sections 
along with a summary discussion. 


C. Radiography. 


I 

1. Introduction, Radiography, as one of the tech-1 
niques used in the satellite inspcaion program, can give a \ 
large amount of infornution concerning size, shape, and j 
density distributioiL of internal components of an inter- 
rogated satellite j 




The basic approach in the radiographic program was 1 
to design and fabricate an imaging system approximating I 
a flyablc system that could easily be used in a scries of j I 
experimental parameter studies. The parameter studies jl 
would define a minimum energy and intensity that would |l 
produce interpretable radiographic images of nuclear ^ 
weapons, as well as investigate target-associated scattering, 11 
the effects of spacecraft structure and of shielding coun- I ^ 
termeasures, and the optimization of components. 

The system that was designed and assembled is \ 
described in some detail in Appendix B, along with the 
results of component evaluation studies. With this system 
a variety of nuclear weapons were radiographed under 
widely varying conditions. Radiographs of a W-59 were 
taken at LASL using the 2.3-McV Febetron and SCAPP at 
3 to 7 McV. At the SCAPP, the W-59 was also radiograph¬ 
ed in a variety of orientations, and sheets of absorber 
material were introduced in cront of the weapon to simu¬ 
late spacecraft structures and/or countermeasures. Experi¬ 
ments using radioisotopic sources were also conducted 
with reasonably good results. The system was also used to 
radiograph a Mk-45, Mk-5 3, W-62, and W-68 with the 
Pantex linac, variable in energy from 7 to 10 McV, as the 
radiation source. In all cases, entire weapons were radio¬ 
graphed, although the pictures were in several segments 
when the size of the weapon dictated. Finally, conven¬ 
tional x-ray film radiographs were taken at all of the 
above radiation facilities and at the LASL betatron 
(22 McV) and million-volt machine to provide a standard 
for comparison. 


3. Experiments at the Febetron. The Febetron is a 
2.3-McV flash x-ray machine that has a nearly 


*Dose measurements throughout this study were made with 
thermoluminescent detectors (TLD's). 
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Scan of W-59 using ^Co. 


monocncrgctic electron beam with a pulse width of 
^ 20 nsec. After several radiographs were taken with both 
the image amplified system and with film, it was 
determined that the machine emined a large amount of 
radiation not associated with the beam spot. A conical 
collimator, together with some lead shielding around the 
face of the machine, reduced this background to a 
tolerable level. Throughout the experiments a 


0.025-in.-thick tantalum target was used directly behind 
the conical collinutor and the beam spot size was defined 
at 1/4-tn. diam. 

Some beam filtering was tried in an attempt to 
harden the spectrum. Slightly better penetration was 
observed with a 1/8-in. lead filter at a cost of roughly 40% 
of the beam intensity. Several of the radiographs taken 
are shown in Fig. 20. For this sequence of shots, the 
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Fig. 20. 

Scan of W-59 using the Febetron. 2.3 MeV. 



sourcc-co-targct distance (Di) was 22 ft and the 
targct-to-scrccn distance (D 2 ) was 5 ft. The dose on the 
screen was 3 mR. The Drincinal details of the device 


t 


^are all visible in these shots. 

^ A senes ot shots was made with various thicknesses 
of aluminum plate which simulated spacecraft structure. 

The shipping container was also present in all of the W-59 | 

radiographs. Figure 21 shows a sequence covering^ T / 

y fhcTirst picture has no ihtcrvcmhgTttenuatoF^^ 

^he second has a l/4-in.-thick sheet of aluminum, the 




third has a l/2-in.-thick sheet, and the founh has a 
l-in.-thick sheet. The 4- by 6-ft sheets were directly in 
front of the W-59. The effects were minimal even at this 
low energy. _ 
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Fig. 21, 

Simulated spacecraft structure, Febetron, 2.3 Mev. I 


4, Experiments at the Single Cavity PHERMEX 
Prototype (SCAPP). The SCAPP is a flash x-ray source 
variable in energy from 3 to 7 MeV with a 0.2-Msec pulse 
width and an almost monochromatic electron beam.* 
The spot size on target is controlled by focusing lenses 
near the target and by beryllium internal collimators, 
thereby allowing the spot size to be varied to suit the 
experiment. Of course, as the spot size is decreased, a 
penalty is paid in beam intensity on target. The 
radiographs taken with the SCAPP were with a 3-mm spot 
and with the conical collimator described earlier. Only 
two targets were used for all the experiments: a 0.010-in. 
tungsten target at 3 MeV and a 0.043-in. tungsten target 
at 4 to 7 MeV. 

Throughout the experiments at the SCAPP, the 
image intensifier was operated either in a dc mode, where 
the tube was on for seconds; or in a pulsed mode, where 
the tube was on for the duration of a gate pulse 
synchronized with the beam burst. Several radiographs 
were t aken using from ^2 tq ^i 2 _ 2 ulscs^ 


Ml 


000 ^ 


•The average current per pulse on target at 3. 4, and 7 MeV 
is 2.8, 6.6, and 21.6 A, which corresponds to 1.68, 5.28, and 
30.2 J/pulse, respectively. 
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intensity on the screen varied from 1 mR (Dj + D 2 =7 ft) 
to 60/iR (Di + D 2 *28 ft). This sequence of photos is 
shown in Fig. 23(a-e). Figure'23(a) shows a view of the 
primary end with a dose of 1 mR on the screen. Figures 
23(b) and 23(c) arc two differ ent views with 25 0 /iR . 
with 23(c) showing the] 


’ I _^_r— 

I Figures 23(d) and 23(c) show the same region with 100 

! and 60 /iR on thc screen rcspeaivcly. As can be seen] 
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Setails visible inside the radi atipn 


4 MeV. A scan at 4 McV, taken with a 0.043-in. 
tungsten target,, and with Dj =» 10 ft and Di -*• D 2 = 12 ft, 
is shown in Fig. 24(a-d). The dose on the screen for this 
sequence was 3 mR.j 




m. 


. i? 




3 McV. A scan of the W-59 taken at 3 McV, using a 
0.010-in. tungsten target, with Dj * 5 ft, = 2 ft and 
3 inR on the screen, is sl^wn in Fig. 22(a-d).| 
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A scries of radiographs was taken at 3 McV with the 
sourcc-to-target distance (Di) varying from 5 to 26 ft. 
^ Because the beam intensity remained unchanged, the 
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Fig. 22. 

(a-d) W-59, 3 MeV, (e) W-59, 3 MeV, 1-in. aluminum, simulated spacecraft structure. 
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It is interesting to note that the Mk-5 3, ] 

radiographed at Pantex at 7 and 10 MeV, weighs almost 
tSOOO Ib) 
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A series of radiographs was made in which the 
intensity on the screen was varied from 1 mR to 15 fJtR, 
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Fig. 24. 

ScanofW-59, 4MeV. 


si milar to the series made at 3 McV.^ _ _^ 
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The effect o? source-associated bac%rouncl sniluld ' 
decrease with decreasing source-to-screen distance. Figure 
27(a-e) illustrates this point. This sequence of pictures 
was taken wit h Di « 3-1/2 ft, D 2 * 1-1/2 ft, and 10 mR 
on the screen.jT V 




5 MeV. Several scans of the W-59 were taken at 
5 MeV using the 0.043-in. tungsten target, with 
Di = 26 ft and Di -»■ Dj * 29 ft. T he dose on the scr een 
^iot Tffgy y mP j ^ ^ 


kO'j 



7 MeV. The scan shown in Fig. 28(a-d) was taken 
with the 0.043-in. tungsten target at 7 MeV. with 
Di = 26 ft, » 3 ft, and with^thc dose o n the sgeen 
equal to 5 mR.J_' _ __ 
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Fig. 25. 

Scan ofW-59, 4MeV, 1.8 magnification. 
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fw simulated ^spacccratt^and b'iglit shing 
experiments were done with Di = 26 ft and = 3 ft» 
similar to those d escribed above using 1 -in . aluminum and 
Hn. steel sheets J 


} 


\ 


^utulTCT^ 


bkcbassitied 


hO) 


31 














































i. 

( 

i 

\ 

I 

I 

j 


I 


Fig. 26. 

W-5P, attenuation with aluminum and steel, 4MeV, 
(a) No absorber, (b) I-in. aluminum, (c) I-in. alumi¬ 
num plus I -in. steel. 


An unattenuatcd sequence was also taken with 
Di = 26 ft and D 2 = 39 ft. The W-59 was scanned in six 
frames as shown in Fig. 29(a-f). If this sequence is 
compared with the sequence shown in Fig. 28, it is 
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Fig. 27. 


Scan of W-59 showing effects of reducing source-associated background, 4-MeV. Di = 3-1/2 ft, 
Dj = 1-1/2 ft. 
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Fig. 28. 

Scan ofW-59. 7MeV. 
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Fig. 29. 

Scan ofW-59, 7 MeV, 1.3 magnification. 
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Tabic X summarizes the various energies used to 
radiograph the W-59 using the image-amplified system and 
the qualitative results for penetration and definition of 
the principal weapon characteristics. The results of 
varying the intensity on the screen at a fixed energy arc 
similarly summarized in Table XI. 

5. Experiments at the Pantex Linac. The Pantex 
Unac is a variable energy (7- to 10-McV) Uncar accelerator 
with a monochromatic electron beam designed specifi¬ 
cally as a radiographic source. Source-associated back¬ 
ground from the Unac was relatively low, but the room 


scatter contribution was more significant. The linac is 
housed in a concrete building with earth shielding on 
three walls. The maximum attainable sourcc-to-scrccn 
distance was 25 ft, but most of the radiographs were 
taken with a source-to-scrcen distance of 19 ft. 

The image intensifier was operated in both the dc 
and the pulsed mode. Because of the very high dose rate 
even at 15 pps, very short exposure times were used in the 
dc operation of the tube. When the image intensifier was 
used in the pulsed mode, a trigger pulse from the Unac 
was used to turn on the image tube at the appropriate 
time. 
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Fig. 30. 

Vif-59, magnifications at 7 )AeV. 


Mk-53. The principal reason for our going to Pantex 
was to use the image-amplified system to radiograph the 
Mk-5 3. This weapon is the largest currently in stockpile, 
the basic weapon weighing 6750 lb. During these 
experiments, the weapon was mounted on the normal 


shipping truck which consists of a thin steel overcoat and 
a heavy structural cart (8000 lb total). 

Two scans of the device were taken with D| « 14 ft 
and Dj * 5 ft at 7 MevJ 
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Fig. 32. 

Scan ofW-59, oblique angle experiment, 7Mei'. 
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Fig. 33. 

Scan of W-59 showing effects of reduction of the source^assodated background, 7MeV. Dx = 12 ft, 
D^-Sft. 
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Fig. 34. 

Scan of W-59 showing effects of further reduction of source-associated background, 7 MeV, (a-d) 
Di = 5-1/2 ft, D 2 = 1-1/2 ft. (e) D, = 3-1/2 ft, = 1-1/2 ft. 


Mk-45. The Mk-45 is a small single-stage weapon, 
and is included here only as an example. The total weight 
of the device in its shipping container is approximately 
150 Ib. Figures 37(a) and 37(b) are two radiographs taken 
at 7 MeV and 30 mR on the screen but with different 
geometries. Figure 37(a) was taken with Di » 15-1/2 ft, 
D 2 » 3-1/2 ft, ^d^ig. 37(W h^ Di= 15-1/2 ft and 
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W-68. The W-68 is a small two-stage weapon used as 
the nuclear device in a MIRV system on the Poseidon 
missile. Figure 38(a-c) shows three views of the weapon 
taken at 7 MeV and 60 mR on the screen and with 
Dj = 15 ft and » 4 ft. In Fig. 38(d), is 10ft, 
slightly enlarging the device. The bar showing at the 
center of the secondary is a part of the shipping and 
handling container. The outstanding features of this * 
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Source E 

192 „ 
<0.600 
MeV 

'^ 1.25 

McV 

Fcbctron 
2.3 McV , 


3.0 McV » 


SCAPP 
4.0 MeV 

SCAPP 
5.0 McV 


SCAPP 
7.0 McV 


TABLE X 

ENERGY VS PENETRATION 
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TABLE XI 

INTENSITY VS PENETRATION 


Energy 

(McV) 

3.0 


4.0 


7.0 


Intensity 

(mR) 


1 

0.25 

0.100 

0.050 

1 

0.063 

0.015 

1 

0.100 

0.050 


Definition of W- ;9.. 
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Mk-53, 7-MeV scan, 2mR. 



HtT comparing Fig s. 38(c) and 38(d) (D^ = 4 and 
10 ft respectively)^ , 
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^j only two views arc shown in Fig. 39, with 39(a) 6. Radiographs with Standard X-Ray Film. 

taken at 7 MeV and 39(b) at 10 MeV, each with 60 mR radiographs were made with the standard industrial 

on the screen. In both pictures, Di = 15-1/2 ft and techniques commonly used at such facil^ ,, ^ 
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Fig. 37. 

_ Mk-45, 7-MeV scan. 

LASL Group GMX-1. Several sources were used in making 
these pictures, including the million-volt machine 
(1 MeV), a 600-mCi cobalt-60 source (1.25-MeV average 
energy), the Febetron (2.3-MeV flash x-ray machine), the 
SCAPP (3- to 7-MeV flash x-ray machine), the betatron 
(22 MeV), and the Pantex linac (7 to 10 MeV). The 
Mk-5 3, W-62, and W-68 were radiographed at Pantex, and 
a number of W-59 pictures were made at LASL. 

These radiographs are not reproduced here but 
served several functions during the program. They 
indicated the quality that can be achieved by using the 
standard industrial techniques, served as a way of 
evaluating the performance of the image-amplified 
system, and were useful initially in establishing 
experimental conditions for the numerous sources. 
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Fig. 38. 

W-68, 7-MeV scan. 
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Fig. 39. I 

W-62, (a) 7~MeV scan, (b) W-MeV scan. ^ 


7. Radiographic Image Enhancement. During the 
course of the program, the possibility of enhancing 
radiographic images by computer processing was consider¬ 
ed. Some effort was expended in this area and is discussed 
in Appendix F. 

D. Photonuclear Interrogation. 

i. Radioisotope SourcesA 1 



"VCunLl 


Hr 
























1 


oronrT 




46 




cussi^. 




LD 


2. Electron Accelerator Source Measurements, 
a. Introduction. A variety of photonuclcar proc¬ 
esses can be exploited for weapons detection if brems- 
strahlung radiation from an electron accelerator is avail¬ 
able as the interrogating source. In addition to (7»n) 
neutrons from and ^Be, photofission of thorium, 
uranium, and plutonium isotopes can be achieved with 
electrons of > 6-McV energy. The latter reactions then 
make possible measurements of prompt and delayed neu¬ 
trons and gamma rays. If electrons of energies >10 MeV 
are available, photoneutron production and/or isomeric 
state gamma-ray emission can be used to detect such 
common materials as lead, bismuth, tungsten, tantalum, 
gold, silver, copper, iron, and aluminum. In addition, the 
pulsed-beam capability of most common electron acceler¬ 
ators allows a measurement of characteristic decay times 
of weapon systems following photonuclear excitation, 
and this can be used to infer design information. 
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TABLE XII 


Warhead 


(7.n) COUNTING PARAMETERS FOR 500-Ci SOURCE 



Net Counts and 

Net Counts and 

Average Neutron 

Yield 

Statistical Error 

Statistical Error 

Energy 

(Mt) 

with 20 cps BKG 

with 200 cps BKG 

(Mev) 

1 

i 

1880± 64 

1880 ± 88 

1.1 

I 

1880 ±63 

1880 ± 87 

^ 1.0 

I 

4520 ± 97 

4520 ±114 

0.05 

i 

1120± 51 

1120 179 

^0.10 

i 

I 

19601 66 

1960 1 89 

^0.05 

1 

1 

3640 ± 88 

3640 1 106 

^^0.05 
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b. Relative Pbotonuclear Yields. A general feature 
of brcmsstrahlung-induced reactions is a steep rise in yield 
as a function of increasing electron energy above the 
reaction threshold. This phenomenon is illustrated in Fig. 
40 which shows the bremsstrahlung-induccd photo neu¬ 
tron yields for a variety of small samples over the electron 
energy range from 5 to 15 MeV. These measurements 
were made at the LASL electron prototype accelerator 
(EPA) with a slab-type neutron detector. Up to about 
10 MeV the yields for D^O and beryllium dominate. 
Above this energy photofission yields dominate. Below 
about 8 MeV the photoneutron yields of such high-Z 
materials as bismuth and lead are negligible if ^Be, ^D,or 
fissionable isotopes arc present. The photoneutron yields 
of intcrmcdiatc-Z materials such as copper and iron arc 
negligible at all energies of interest. 



4 6 8 iO t2 M 

ELECTRON ENERGY ( M«V) 


Fig. 40. 

I Relative bremsstrah lung-induced photoneutron 
yields for se ve ral ^terial ^^^, ^ 
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Pbotoneutron yield as a function of electron energy ELECTRON ENERGY ( M#V) 

for three classes of weapon secondaries. ^2. 


Slab-ratio data for bremsstrablung^nduced photo- 
neutro n reactio ns. 








































TABLE XIII 



APPROXIMATE AVERAGE PHOTONEUTRON SPECTRAL ENERGIES IN McV 
AT THREE ELECTRON ENERGIES 


Secondary System 


= 5 McV 

E, « 7.5 MeV 

E^ * 10 McV 

1.08 

1.19 

1.51 

1.00 

1.16 

1.50 

0.25 

0.78 

1.51 

0.50 

0.89 

1.50 

0.25 

0.40 

0.94 

0.50 

0.72 

1.12 

0.75 

1.05 

1.45 


a ^Hc proportional counter. Generally speaking, 
photoneutron spectral shapes and endpoints are unique 
for each isotope at each electron energy. Table XIV shows 
photoneutron endpoint energies for several isotopes of 
interest at three electron energies. Neutron-induced fis¬ 
sion neutron spectra can be considered the same as photo¬ 
fission neutrons from ^U, or other photofission- 

able nuclei. 

There is a large endpoint difference between and 
’Bc at all electron energies. Considering the problem of 
competing reactions, however, an electron energy around 
5 MeV is most interesting for beryllium-deuterium dis¬ 
crimination. Figure 43 illustrates the endpoint technique 
by showing measured beryllium and small sample D^O 
photoneutron spectra for 5.35-McV electron bremsstrah- 
lung. These spectra were obtained by unfolding ^He pro¬ 
portional counter spectral data measured at the EPA. 
Appendix D gives a description of the unfolding tech¬ 
nique. The two spectra differ considerably in both end¬ 
point energy (1.56 and 3.27 MeV) and shape. The spectral 
shape differences are due prinurily to differences in the 
shapes of the D( 7 ,n) and Be( 7 ,n) cross seaions below 
5 MeV. Althoi^h the D( 7 ,n) cross section in this energy 
range is well known,a literature search indicated that 
the Be( 7 ,n) cross section is not well known above 3-MeV 
photon energy. Competing beryllium photo nuclear reac¬ 
tions such as (1) ’Bc+ 7 -^*Bc**t*n (2.9-McV excited 
state of *Be*), (2) ’Be + 7 -^ ^He + ’ He - 2.46 MeV 
(* He subsequently breaks up spontaneously into ^ He + n 
with a 10“^^ sec half-life), and (3) ’Be 7 ^He 
+ ^He n — 1.57 MeV could also occur and would 


make the prediction of neutron spectra uncertain even 
if the total Be( 7 ,n) cross section were known. 

Figure 44 shows (1) the ^^^Cf spontaneous fission 
neutron spectrum, ( 2 ) a ^’Pu multiplying assembly neu¬ 
tron-induced fission neutron spectrum, and ( 3 ) the 
9.8-McV electron bremsstrahlung-induced photo¬ 

fission neutron spectrum measured at EPA. All three 
spectra were measured with ^He proportional counter 
spectrometers and unfolded using the technique described 
in Appendix D. All three spectra arc quite similar, 
except for differences at low neutron energies (due to 
moderation effects in the multiplying assembly and a 
competing photoneutron reaction in the case of ^U). 


TABLE XIV 

BREMSSTRAHLUNG-INDUCED PHOTONEUTRON 
ENDPOINT ENERGIES IN McV 


Isotope 

Eg = 5 MeV 

E^ = 7.5 McV 

Eg » 10 McV 


1.39 

2.64 

3.89 

’Be 

2.96 

5.18 

7.40 

Pb 

(a) 

0.76 

3.26 

Bi 

(a) 

0.04 

2.54 

238y 

(a) 

> 10 

> 10 

235 u 

(a) 

> 10 

> 10 


^*^Bclow chreshold for reaction. 
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Fig. 43. 

Deuterium and beryllium pbotoneutron spectra, 
S.JS’MeV bremsstrabbing. 


' ^ principal advantage of ^^rtional-countcr 

spectrometers are low gamma sensitivity (successful spec¬ 
tra have been taken in fields of > 10 R/h) and high 
intrinsic neutron detection efficiency 13% at 1 MeV 
for a ^He type is realistic). Several experimental compari¬ 
sons were made between ^ He and CH 4 tubes of identical 
scanering pressures (2 atm CH4 and 8 atm ^He) under a 
variety of conditions. This comparison between ^ He and 
CH4 for monoenergedc 1-MeV neutrons is shown in 
Fig. 45. The ^He tube is approximately twice as efficient 
as CH4 in the maximum recoil energy region, and this 
translates into a correspondingly increased efficiency in 
applying the endpoint technique. The ^He counter also 
has a considerable advantage in gamma sensitivity. In a 
> 5-R/h gamma field, the ^He spectroir^ter could be 
used to take good pulse-height spectra down to about 
0.2 MeV, whereas the CH4 counter suffered appreciable 
pulse-height degradation at 0.75 MeV. Identical electron¬ 
ics were used for these measurements, with high voltages 



Fig. 44. 

Fission neutron spectra for ^U, ^Pu, and ^Cf. 


for the two tubes set at their respective optimum resolu¬ 
tion points. The ^He elastic scattering proportional- 
counter spectrometer was thus selected for a further 
study of the endpoint technique. In addition, the proper¬ 
ties of various types and configurations of ^ He spectrom¬ 
eters were studied analytically using LASL Monte Carlo 
computer codes. A description of the analytical work is 
included in Appendix E. 

The basic deuterium-beryllium discrimination using 
the endpoint technique is illustrated in Figs. 46 and 47. 

Figure 46 shows the measured small sample photoneutron 
^He spectrometer pulse-height distribution for D^O and 
beryllium for bremsstrahlung from 4.5-MeV electrons. 

Figure 47 shows the same for 5.35-MeV electrons. The 
corresponding endpoints are marked in all cases on the 
horizontal axes. Although these are direct spectrometer 
data (no unfolding), the shapes of the deuterium and 
beryllium distributions at a given electron energy arc 
characteristically different, as arc the endpoints. In an 
unknown spectrum with only D( 7 ,n) and Bc( 7 ,n) compo¬ 
nents (c.g., conventional and dean weapo ns i nterrogated 
at less than 5 MeV), the beryllium portion could be iden¬ 
tified and the deuterium component extracted for spectra 
containing a beryllium component less than a few times 
the deuterium component. The beryllium and deuterium 
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Fig. 45. 

^He and CH 4 IrMeV neutron response, neutron 
spectrometer comparison. 



Fig. 46. 

^He neutron spectrometer data, d.S-MeV 
bremsstrahiung-induced pbotoneutron pulse-height 
distribution for beryllium and D 2 O. 



Fig. 47. 

^He neutron spectrometer data, S.SS-MeV 
bremsstrablung-induced pbotoneutron pulse-height 
distribution for beryllium and D^O. 


pulse-height distributions are shown for both 4.5- and 
5.35-MeV electron energies to illustrate the approximate 
range over which a successful discrimination can be made. 
As the electron energy is raised, all yields are correspond¬ 
ingly raised, with roughly a factor of two increase in yield 
per electron between 4.5 and 5.35 MeV for both deute¬ 
rium and beryllium (Fig. 40). However, the ^He detector 
efficiency fails off monotonically for neutrons >1.1 
MeV. This factor must also be considered in arriving at an 
opti mum electron energy. 
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When higher energies are used, the photofission 
neutrons of ^*U, and ^^Pu can be used for weapons 
detection. Figure 49 shows the ^He spectrometer meas¬ 
urements of small sample photofission neutrons. 

Pulse-height distributions are shown for 6-, 7-, 10-, 12-, 
and 15-MeV electrons and are normalized to the same 
approxinute bremsstrahlung dose incident on the sample 
at each electron energy. If the curves are multiplied by 
(Eg)^ at each elearon energy, they will be expressed 
approximately on a relative scale per electron. 

The D(7,n) distribution for 6-to 15-MeV electrons, 
also plotted for equal doses at each electron energy, is 
shown in Fig. 50. The systenutic change of the spectral 
endpoint can be easily seen. Figure 51 shows a compari¬ 
son of and Pb<7,n) signatures at 10 MeV and 



Fig. 48. 

^He neutron spectrometer data, 5.5-MeV 
bremsstrahlung-induced photoneutron pulse-height 
distribution for D^O and for D^O plus assem¬ 
bly. 



Fig. 49. 

^He neutron spectrometer data, photoneutron 
and photofission neutron pulse-height distribution 
for various bremsstrahlung energies. 

illustrates a possible means of distinguishing between con- 
ventional and clean secondary designs^ _ 
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Fig. 50. 

D(y,n) ^He proportional counter pulse-height distri¬ 
bution for 6- to 15-MeV bremsstrahlung. 



d. Calculated and Scaled Results. To some extent 
the various results presented above can be summarized by 
scaling to the same systems described in Table XII for the 
radioisotope (7,n) results. Anticipating the system defini¬ 
tions presented in detail in Sec. VI, two specific source 
and detector combinations are selected. 

Table XV shows the number of counts obtained for 
various weapons using a Marx generator flash x-ray source 
at several energies and the slab deteaor described earlier. 
The Marx generator and the neutron detector arc assumed 
to be 20 ft from the weapon. The net counts shown arc 
for 40 pulses of 1.25 x 10^^ electrons each (each pulse is 
thus equivalent to 100 A for 0.2 Msec). These calculated 
values arc based on small sample measurements, radio¬ 
active source (7,n) measurements, and the attenuation 
and multiplication calculations mentioned above. The 
expected count from beryllium in small RV quantities is 
also shown. For beryllium in this quantity or less, a few 
pulses at 2.2 MeV would suffice to determine beryllium 
effects to ± 20%. However, this technique is not a good 



CHANNEL NUMBER 


12 3 4 

APPROXIMATE NEUTRON ENERGY ( MtV ) 

Fig. 51. 

10-MeV bremsstrahlung^nduced photoneutron 
pulse-height distribution for lead and^U. 


one for larger quantities of beryllium) 

1 

Because Marx generator measurements would be 
pulsed, the effects of space backgrounds would be small. 
Assuming the backgrounds previously discussed and a 
50-/iscc counting time per pulse, the total background 
count from all sources during the 40-pulsc irradiation 
(Table XV) would be less than one count. 

Table XVI shows brcmsstrahlung-induced count 
rates expected for a photointerrogation system using a 
500-W, 5-McV electron linac. This beam is incident on a 
thick (*^ 0.10 cm) tungsten bremsstrahlung target assum¬ 
ed to be 10 ft from the interrogated weapon. The linac is 
farther away with the bremsstrahlung target suspended on 
a lightweight boom. The detector is assumed to be a bank 
of 10 16-atm ^He proportional counters (2 in. diam, 
24 in. long) located 20 ft from the interrogated weapon. 

The entries in Table XVI were obtained by scaling 
measured D(7,n) ^He spectrometer count rates from the 
EPA small sample measurements to the ^ LiD mass distri¬ 
bution appropriate to each weapon. Complicating effects 
such as bremsstrahlung attenuation in intervening high-Z 
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Fig. 52 ' 

Linac interrogation of the Mk~53 for electron ener¬ 
gies of 6 to 12MeV. 

materials and (n,f) reaction components were taken 
into account using the techniques described in Ref. 9. The 
resuiting output spectrum (analyzed approximately into 
separate D(7,n) and (n,f) components) was then 

folded into the ^He detector response, and the final 
spectrometer response was summed over the indicated 
limits to illustrate approximate statistical accuracies for 


various pulse-height ranges. The beryllium-only entry was 
calculated in a similar fashion assuming that 2.5 x 10^ 
cm^ .of beryllium, 0.97 cm thick, is on the exterior por¬ 
tion of the target weapon. 

The effect of space background on these measure¬ 
ments IS negligible, primarily because the linac is pulsed, 
and the indicated 100-sec interrogation corresponds to a 
detector on-time of about 0.5 sec. Using the space neu¬ 
tron background estimates already discussed, about 4 to 
40 total counts can be expected, distributed over the 
entire pulse-height range 0 to > 10 MeV. This background 
could probably be neglected entirely, however, an in situ 
background measurement between beam bursts could 
easily be made. 

liicctronics for the ^He spcctrornctcr arc identical 
to chose discussed earlier for cac ^ He slab dctcaor with 
the exception that a pulse-heignc analyzer is also required. 
Only a modest number of channels of pulsc-hcight data 
arc required 100), and the puise-height analyzer dis¬ 
cussed ui c-onncction with the passive ganuna detector is 
more than adequ^e for this purpose. 



The spectra of Fig. 52 can be converted 
conditions of Table XVI by simple scalingj 




to the / 






TABLE XV 


APPROXIMATE SLAB DETECTOR COUNT RATES FOR 
40 PULSES OF A MARX GENERATOR WITH 
100-A, 0.2-iLiscc PULSES 


Weapon 
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4-McV Electrons 
(Net Counts) 


474 

454 

1140 

282 

494 

918 


474 


5-MeV Electrons 
(Net Counts) 


2.2-McV Electrons 
(Net Counts) 


1452 

1390 

3490 

866 

1514 

2812 


1058 
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TABLE XVI 

APPROXIMATE COUNTS OBTAINED WITH A 100-SEC, 

5-MeV LINAC INTERROGATION OF SELECTED WEAPONS SYSTEMS 


Weapon System 




Counts > 1.39-MeV 

Counts > 2.96-MeV 

Counts > 0.5- 

Counts >1.0- 

Pulse Height 

Pulse Height 

MeV Pulse Height 

MeV Pulse Height 

[D( 7 ,n) Endpoint] 

[Bc(7,n) Endpoint] 

5000 

1970 

960 

120 

7300 

1000 

35 

4 

4800 

1890 

920 

no 

1800 

250 

9 

1 

5880 

800 

28 

3 

1000 

260 

70 

— 
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ments to compare relative 14-MeV neutron and 10-MeV 
electron bremsstrahlung-induced fission efficiencies were 
performed on a Mk-5 3 at Pantex ii^a standard geometry 
\ using the same neutron detector. I 


a- 






4. Delayed Fission Neutron Time History. Delayed 
neutrons are a unique signature for fissionable isotopes at 
10 MeV because there are no other photo reactions on any 
isotope that can generate directly a delayed neutron. In 
addition, under favorable conditions the delayed neutron 
time history from photofission can be used to determine a 
rough isotope ratio. 

Figure 5 3 shows a delayed neutron time history for 
the Mk-5 3 system foUpwing 10-MeV bremsstrahlung- 
induced photo fi ssion, y "" ^ ^ 


-- I 


Figure 54 shows the delayed neutron time history I 
for a 7-MeV bremsstrahlung irradiation of the XW-67 
secondary. This measurement was done at the SCAPP at 
LASL. and the data shown are the sums produced by 10 
O.I-fisec pulses spaced 30 sec apart. Again, a slab detector 
was used. The solid line through the data points is the 
calculated time history using published <n.f) delayed 
neutron time history parameters.The dashed line is the 
calculated time history using published small sample 


I 

I 

IPTg. 53. 


Mk-53, 10-MeV bremsstrahlung-induced delayed 
neutron time history. 
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__ Because the space background can be 

nieasured in some cases, a dme history of sorts can 
appar ently be m easu^d fpr_jaige^ weapo ns at 10 MeV. 




5. Time-of’Fligbt Spectrometry, If an inspection 
could be carried out at distances of 5*100 ft under stable 
ranging conditions (AR/R < 0.05) with short beam bursts 
(< 0.1 /isec), then a time-of-flight technique coidd be^ 
attemp ted/ 


Fig, 54, 

XW-67 secondary, 7-MeV bremsstrabhtng-mduced 
delayed neutron time bistory, 

(y,f) delayed neutron time-history parameters.^^ An inde¬ 
pendent measurement done at SCAPP using a snull U 
sample verified the values in Ref. 21. Figure 54 demon- 
y V strates tlwt even with 7-MeV bremyrah lung. tl^ time 
^ history | 


'd 




'Sinular time histories were taken 
for tfic XW-67 at 4, 5, and 6 MeV, and all displayed the 
shape shown in Fig. 54. 

The delayed neutron measurements of the Mk-5 3, 
XW-67, and small samples of ^®U and lead one to 
expect at least two distinct delayed neutron time histories 
from weapons following interrogation in t he 4- to 10-MeV 
bremsstrahlung energy range.! 


i 


(There is a factor 




of 1.8 difference in an early-to-late ratTo (count rate at 
^0.5 sec compared to count rate at 20 sec) between 
these two delayed neutron time histories for identical 
short burst irradiations. However, it is not likely that 
conditions of a typical space inspection mission will be 
sufficiently stable to exploit this difference. 

The 10-MeV Mk-53 delayed neutron experiment, 
when scaled to a 1-sec irradiation of 5 x 10"^ A of 
lO-MeV electrons incident on a bremsstrahlung target 
10 ft away with the slab detector 20 ft ^om the Mk-5 3^ 
results in a delayed neutron count rate 


E. Neutron Interrogation. Experiments and calcula¬ 
tions were done to investigate the feasibility and useful¬ 
ness of active 14-MeV neutron interrogation in the satel¬ 
lite inspection problem. Selected neutron sources and 
assemblies were used to measure prompt decay constants, 
integrated prompt responses, and delayed neutron yeields. 

I, Experimental Results. Measurements were made 
at Pantex on the Mk-53 and the W-68 using a Sandia 
prototype, high-output neutron generator. The weapons 
were pulsed once every 300 sec for a total of three pulses 
of approximately 10^® neutrons (14 MeV) per pulse to 
provide both delayed neutron kinetic response and delay¬ 
ed neutron yield data. 

The detector consisted of a slab of polyethylene 3 
by 20 by 24 in. with 19 1-in. diam ^He proportional 
counters. This particular detector was designed for labora¬ 
tory use but does permit extrapolation to a suitable 
detector for use in space. Pulses from the detector were 
fed through an appropriate preamplifier, amplifier, and 
discriminator. The resultant logic pulses were then fed 
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into a multichannel analyzer operating in the muiriscaic TABLt XVlil 

mode. A crystal-controlled external timer was used co 

provide varying channel widths of 0.5, 1.0, and 23 see EXTRAPOLATED PANTEX OAlYi TOR 

A SOURCE OF 10 PULSES PER SECOND 
OF 5 by lO’ NEUTRONS PER PULSE 

Count Kate 



TABLE XVII 


PANTEX DELAYED NEUTRON YIELD DATA USING 
SANDIA PROTOTYPE HIGH-OUTPUT NEUTRON GENERATOR 


Weapon System 


Distance to Center 
of Weapon 
(in.) 


Interrogation 

Location 


Net Count Rate 
(from t = 0 
to 0.5 sec) 


Source Detector 
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Mk-53 delayed neutron decay following 14-MeV 
neutron interrogation. Sandia h^b-output neutron 
generator. 


2. Calculated Decay Constants and Multiplication 
Factors. The decay constant, a, and the effective neutron 
multiplication factor, were calculated for all weapon 
prinuries listed in Table XIX. These calculations were 
made with the LASL DTF-IV neutron transport 
program^*^ in spherical geometry. The Sg discrete ordi¬ 
nates approximation was used with 16 energy groups and 
isotropic scattering. Cross-section values from the current 
LASL 16-group Hansen^^ set were used. 


Calculation of for each of the subcritical prima¬ 
ries is straightforward and involves no unusual computa¬ 
tional problems. Calculation of the decay constant, a, 
involves the basic physical assumption that the time be¬ 
havior of the spatial and energy dependent neutron den¬ 
sity, n(r,n,v,t), can be described by a fundamenui decay 
constant, oe, such that 


n(r,v,n,t) = n(r,v,12)exp(at) , 


( 1 ) 


where a < O for the subcritical systems considered here. 
The assumption of a dominant fundamental decay con¬ 
stant implies that higher order spatial modes and kinetic 
distortion of the neutron density will decay much faster 
than the fundamental mode and the refore will not, 
to the m ea sured decay rates. [ 

I 


The method of computing decay constants will be 
discussed briefly. Using standard notation, the neutron 
transpon equation is written as 

^-n(r,v,n,t) 2 (v)0(r,v,n,t) r2'V0(r,v,n,t) 

ot c 


= S(r’v,f5,t) 


f 0(F,v\Jl\t) 


2 (v')f(v,r2; v',n')dv' dn' 


( 2 ) 


^-n(r,v,n,t) = ^0(r,v,f2,t) exp(at) 


(3) 


y 


The transport equation is reduced by the method of 
discrete ordinates^^ to a set of coupled multigroup differ¬ 
ence equations. These equations are ordinarily solved 
numerically for the cigejivaiuc, k^f^, which is contained in 
the source term S(r,v,^2,t). To calculate the decay con¬ 
stant, a, Eq. (1) is substituted in the first term of Eq. (2), 
resulting in 


Because the time dependence exp(ort) is a multiplier in 
each term of Eq. (2), the transport equation reduces to 
the steady state expression 

-J- f2-70(r,v,n) = S(r,v,r2) 
0(Fv',n') 25 (v')f (V.H; v'.n*) dv' . (4) 


Equation (4) can now be reduced by discrete ordinates 
techniques and solved numerically with a treated as the 
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__ These instabilities are reasonably^ 

‘ well understood and can be overcome by repeatedly 
adjusting convergence rate control Darameters. Calculated 
decay constants in Table XIX are definitely converged 
values and should contain numerical uncertainties no 
larger than a few percent. /\greement between calculated 
and measured a values is generally good and indnates that 
the techniques used to measure deca y constants are accu- 
rarg* ixnr\ fiMMlCn** j 


CaTcuiations’^' lifdicate that it mlT be very difficult 
or impossible to measure both the prompt fundanwptal 
V, n^de and the prompt-to-delayed ratio\^ 

^ —, 


Me OcfciTom 

f .* 


Calculated Value 
Using l*Oi<ncnsionRl 
DTI*-!'. 

(set * 


Effective 
Muitiplicadon 
Factor, k^£| 
Calculated Using 
DTF-IV 


\ 


I 


i. Neutron Source. Sandia Laboratories (Albuquer¬ 
que) are developing a high-output, 14-MeV neutron gener¬ 
ator lO produce liie order of lO^n/Msec in pulses with 
widths of 5 to 30 Msec. Presently they have atuined 
3xlO*n/Msec with generator lifetimes of greater than 
400 pulses. The high voltage and ion source are on only 
during the pulse, because a pulse transformer is used to 
generate the source and acceleration voltage. Thus, the 
background between pulses is only the result of activation 
of in the tube insulators, and these background levels 
are very low. For the application in this study, i.e., the 
measurement of prompt responses of weapons, the 5-Msec 
pulse would be used. 
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The neutron generator and pulse transformer can be 
packaged in a cylinder about 15 in. long and 7 in. in diam 
and weighing about 15 lb. The power required is approxi- 
nutely 200 mA at 4 kV for the duration of each 20-sec 
active interrogation. In addition, a few wans (less than 
50) would be required during the active interrogation and 
the background determinations. 

The deteaors chosen were a cluster of five proton 
recoil detectors with minimal shielding and the small slab 
detector as described in the (7,n) technique. The inspec¬ 
tion distances are not great enough to cause distortion of 
the time spectra for neutrons above 0.5 MeV from differ¬ 
ences in dme-of-flight. No gamma-ray discrimination is 
used because it is difficult to operate at the high count 
rates, and a higher efficiency is obtained by counting both 
neutrons and gamma rays. The suitability of this mode of 
operation can be seen by observing the agreement for d ie 
decay constants given | ^as 

shown in Table XVIII. ^he proton recoil detectore were 
chosen over ^ He or CH^ proportional counters because of 
the increased efficiency attainable. The small slab de¬ 
tector used in the (7,n) technique will be used here for 
the detection of delayed neutrons. 

The devices in Table XX give a representative sam¬ 
pling of the various types of weapons that might be 
encountered and thus represe nt reason^le limits on 
count rates. If a system | 


TABLE XX 

ESTIMATED INTEGRATED COUNT AND THE 
ASSOCIATED EARLY AVERAGE COUNT RATE 


System 


/c 


JCdt 

20 fJLSCC 


^T^ 20 Msec 

JCdt 

20 Msec_ 

■■vi 




3 




the detected count rate will be re duced 
considerably. The sa me is true for \ 

4 In thi^ase the 

systems may not bTiietcctablc witli this method. 

F. Delayed Gamma Rays. An inspection system 
which has a pulsed neutron or photon source on board, as 
well as a Ge(Li) detector, can obtain additional infor¬ 
mation by analyzing the gamma-ray spectrum between 
interrogating pulses. There are several potential sources 
for the gamma rays emitted following an incident pulse of 
photons or neutrons. 

Delayed fission gamma rays are emitted by fission 
products which are left in high angular momentum states 
after fission. Lifetimes are long 10“^ sec) because tran¬ 
sitions from these high angular momentum states are 
highly inhibited. For a weapon system, experiments indi¬ 
cate that these delayed fission gamma rays are extremely 
weak relative to gamma rays from the other sources to be 
discussed. 

Isomeric states in nonfissionable materials can be 
excited via (7,n), (n,2n), and possibly (n,n') reactions. 
Over the periodic table, the lifetimes of these states vary 
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for a 3-in .-chick mock HE shell with an incident pulse of 
10^ fast neutrons, the thermal neutron flux in the shell 
was '^10^ neutrons/cm^/sec and had a half-life of 
1 msec. This flux estimate is based on measurements of 
the 2.224-MeV gamma ray from the H(n,7)D reaction in 
the shell. During this thermalization and die-away, (n,7) 
and (n,f) reactions can take place. The neutron capture 
will produce characteristic gamma rays observable with a 
Ge(Li) detector. These gamma rays will have energy, 
relative intensity, and half-life (which will follow the 
neutron half-life) signatures which should eiuble identifi¬ 
cation of the elements in the hydrogenous material. For 
example, in fast neutron interrogation of the mockup 
systemsj gamma rays from H(n, 7 ), Q(n,7), 

and N(n,7) were clearly identified. Direct evidence for 
fissionable materials was not found except that in some 
cases prompt fission neutrons caused Ge<n,n'7) reactions 
in the Gc(Li) detector. At 690 keV there is a peak which 
has a very characteristic shape. 

The (n,n') and (n,2n) reactions, which will occur 
during the neutron burst, can produce nuclei in isomeric 
states. The characteristic gamma rays from isomeric states 
in “Pb and *^Pb have been observed.! 


from <10“’ sec to tens of minutes. Many of these iso-I 
mers emit gamma rays with characteristic energies, rela-j 
rive intensities, and half-lives. For example, the (7,n) or) 


209t 


208 . 


Bi which has a 2.5-msec half-life. This state emits 
0.925-, 0.654-, and 0 JlO-MeV gamma rays with a strong 
relative intensity and some other gamma rays weakly. The 
same reactions on natural lead yield isomeric states in 
, = 0.13 msec) and = 0. 8 sec), each of 

emiu several gamma rays.[ 

A quick survey of the chart of the nuclides indicates 
that there are 45 nuclei with isomeric levels which have 
half-lives between ICT* and 10 sec, and which can be 
reached by (7,n) and (n,2n) reactions on naturally- 
occurring isotopes. Similarly, there are 45 nuclei which 
have isomeric levels in the same half-life region, and which 
may possibly be excited by inelastic neutron scattering. 

If an interrogation is made with neutrons and there 
is a substantial amount of hydrogenous material present, 
then another process will occur. There will be consider¬ 
able neutron thermalization in the hydrogenous material 
over a period of time like 10“^ to ICT^ sec. For example. 


VI. System Definition 

A. Approach to System Definition. Based on the 
techniques described in the previous sections of this 



















TABLE XXI 

ISOMER COUNT RATE AT 15 MeV 



<1 


Isomers 

Ey (MeV) 

Photopeak Counts 

^Pb 

1.066 

10® , 


0.572 

2 X 10® j 

“•pb 

0.925 

2 X 10® . 


0.510 

3 X 10® 1 


report, several sensor package concepts have been chosen 
as reasonable candidates for a spacebome satellite inspec- i j 
tor. These represent somewhat arbitrary but representa- } i 
tive steps in the continuum of options that develop I 
around the possible variations in weapon design, back- | j 
ground conditions, countermeasures, and detection confi- ! 
dencc. j | 

The following assumptions about the overall inspec- j j 
tion system influenced the conceptual design of the basic | i 
nuclear sensor packages: 
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TABLE XXII 

SUMMARY OF SENSOR PACKAGES 
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Schematic of gamma spectrometer. 
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Fig. 59. 

Schematic of radiographic imaging unit. 
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TABLE XXVI 

DETECTION CAPABILITY OF RADIOGRAPHIC SENSORS 
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INSPECTOH 



(0) SLAB NEUTRON DETECTOR TOP VIEW 
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Fig. 60. 

Source-target’detectOT geometry for radiographic 
and (y,n) techniques. 



(b) SCHEMATIC OP COUNTING ELECTRONICS 



HV 


Fig. 61. 

Slab neutron detector and a schematic of counting 
electronics. 
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TABLE XXVll 

WEIGHT (ESTIMATED) SUMMARY OF 
PROPOSED SENSOR PACKAGES 
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TABLE XXVIII 

DATA HANDLING SUMMARY 


Radiographic sensor 
(vidicon)^*^ 

Gamma-ray spectrometer 
< 2 

He neutron detector 

i 

j ^ He neutron detector 

j Proton recoil detector 
\ 


Information 

Content 

Per Record 
(bits) 

Maximum 

Transmission 

Rate to DHS 
(bps) 

8.25 X 10* 

24.75 X 10* 

1.60 X 10* 

8 X 10* 

1.7 X 10* 

1.7 X 10* 

1.9 X 10* 

9.5 X 10* 

5.2 X 10* 

3.6 X 10* 
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APPENDIX A 

Gc(U) DETECTOR CHARACTERISTICS 


Very often a computed or measured count rate is 
quoted for a Ge(Li) detector whose only given specifica¬ 
tion is active volume, i.e., 50 cc. 80 cc, etc. With only this 
specification one cannot, in general, extrapolate the count 
rate accurately to other detectors or geometries. Addi¬ 
tional data are needed concerning the exact dimensions, 
drift depth, and geometry because no two detectors have 
the same geometrical parameters. 

The detector industry now standardizes the specifica¬ 
tion of Ge(Li) detectors by quoting efficiency relative to 
that of a 3- by 3-tn. Nal detector at a source-detector 
distance of 25 cm and an energy of 1.33 MeV. The 
1.33-MeV efficiency for a 3- by 3-in. Nal detector under 
these conditions is 1.2x10"^ photopeak counts per 
source disintegration, a calculated number. A somewhat 
different number would be measured from a ^Co source 
due to summing effects. This convention is followed in 
subsequent discussion; thus the efficiency of a Gc(Li) 
detector expressed as a percent, i.e., 6%. 10%, etc., is 
relative to a 3- by 3-in. Nal detector at 25 cm and 
1.33 MeV. This relative efficiency would be different if 
measured at other source energies or other distances be¬ 
cause of geometrical and photon cross section differences 
between germanium and Nal. The largest Ge(Li) detectors 
now available off the shelf have relative efficiencies of 
about 20%. 

The detector geometry is also important. For the 
purposes of this study a closed end coaxial geometry is 


preferred, because the entire frontal area is aaive. The 
ratio of active-to-dead volume should be as large as 
possible and is well over 95% for state-of-the-art de¬ 
tectors. The drift depth and length-to-diameter ratios are 
important for our application. Larger drift depths increase 
the peak-to-Compton ratio as well as the efficiency per 
unit volume. Drift depths approaching 20 mm arc avail¬ 
able at the present time. 

The interplay of these factors is illustrated in 
Table A-I. Detector 1 is the 6% Ge(Li) detector used to 
record the spectra displayed previously. Detector 2 is a 
large volume state-of-the-art detector of the type pro¬ 
posed for this mission. The volume of detector 3 is even 
larger than that of detector 2 but detector 3 has a much 
poorer efficiency. The ublc vividly illustrates the pitfalls 
of specifying a detector only by volume. The superiority 
of detector 2 over detector 3 is due to greater drift depth 
and a larger diametcr-to-length ratio. Because of its 
length, detector 3 might exhibit better performance for 
this application if used in side incidence. 

Measured efficiency curves for detectors 1 and 2 are 
shown in Figs. A-1 and A-2. These efficiencies were meas¬ 
ured for a source-detector distance of 1 m, the largest 
distance practical with the calibrated source intensities 
available. The large distance was used to approximate the 
nearly normal incidence of photons on the detector face 
that would be encountered in an inspection. The efficien¬ 
cies for a 3- by 3-in. Nal detector shown in these figures 
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TABLE A-I 

CHARACTERISTICS OF SEVERAL Gc(U) DETECTORS 


Detector No. 

Active Volume 
(cm^) 

Diameter 

(mm) 

Length 

(mm) 

1 

35 

40.5 

28.5 

2 

85 

47.8 

48.9 

3 

88 

41 

71 


Drift Depth 
(mm) 

25-cm, 1.33- 
MeV Relative 
Efficiency 
(%) 

Inferred 25-cm, 
1.33-MeV Intrinsic 
Photopeak 
Efficiency 
(%) 

16.5 

6 

4.4 

18.5 

18.2 

9.5 

15 

10.2 

7.3 



Fig.A-L 

Measured photopeak efficiency for gamma spec- 
tro meters. 



Fig. A-2. 

Comparison of gamma spectrometers. 


are calculated and measured values taken from Phillips 
Petroleum Co., Idaho Falls, report IDO-16880 (August 
1964). In these figures intrinsic efficiency is defined as 
photo peak counts per gamma incident on the detector 
face. From the curves of Fig. A-2, which have the de¬ 


tector area folded in, count rates can be computed 
directly when the flux is given. Indeed, the data in 
Fig. A-2 were used to compute the flux from measured 
photo peak count rates, and also to predict count rates for 
a large detector when the flux is given. 
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APPENDIX B 

THE RADIOGRAPHIC SYSTEM 


1. Description 

The basic radiographic system is a 
phocofluorographic system using a scintillating screen 
coupled by an objective lens to an image intcnsificr. The 
image output of the intensificr is then recorded directly 
on Polaroid film through a fiber-optics coupling plate. A 
3- by 3-ft scintillating screen was selected as being conven¬ 
iently large for the radiographic targets of interest and 
sufficiently portable for transporting to a variety of radia¬ 
tion facilities. The screen composite is made up of 
200 mg/cm^ ZnCdS (US Radium PFG) covered by a sheet 
of lead 0.060 in. thick and a sheet of aluminum 0.060 in. 
thick. A lighttight bellows construction was selected to 
permit the use of a variety of lenses with different view 
angles. Most of the experimental work was done with 
only two lenses, although five lenses were evaluated and 
will be described later. 

The image intensifier is a second generation prox¬ 
imity-focused intensifier which is very simple in design 
(although complex from the manufacturing point of 
view), very small and lightweight, and which requires 
relatively low supply voltages (about 6000 V) and very 
little power (about 6 mW). Two intensifiers were pur¬ 
chased from ITT for this program, one of 18-mm format 
and one of 23-mm format. The specifications for the 
tubes are given in Table B-I. 

The 18-mm tube is of special interest because it can 
be coupled to a standard size vidicon target (also 18 mm). 
For imaging directly on film, the larger, 2 3-mm tube was 
more desirable. Second generation tubes as large as 
75 mm have been produced, but such tubes are currently 
produced in small numbers with long delivery times. 

These image intensifiers have a noise-induced out¬ 
put, with no light input, much like photomultipliers. 
Unlike photomultipliers however, there is no convenient 
way to discriminate against the noise. Therefore, when 
the image is being recorded directly on film, the noise- 
induced light is being integrated by the film in just the 
same way as the radiographic image. This noise then 
imposes a time limitation on the exposure time when 
isotopic sources are being used. It can be partially or 
totally eliminated by '*gating’* the tube on and off when 
using a pulsed source, providing the noise is produced in 
the photocathode. 

A gating system was designed to switch the image 
intensifier rapidly, since the tube takes approximately 
1 sec to turn on after full voltage is ap'plied. The tube 


TABLE B-I 

IMAGE INTENSIFIER SPECIFICATIONS 



18-mm Tube 

23-mm Tube 

Maximum usable 
aperture 

18 mm 

23 mm 

Maximum diameter 

45 mm 

50 mm 

Thickness 

22 mm 

17 mm 

Conductivity 

(amperes 

at minimum gain) 

6.7 X 10"“^ 

1.6 X 10’* 

Phosphor 

P-20 

P-20 

Photocathodc 

S-25 

S-20 

Resolution 

24 line pairs/mm 

20 line pairs/mm 

Cathode sensitivity 

232 /iA/lm 
at 4.3 X 10“^ 

Im 

293 *iA/lm 
at 7.5 X 10'^ 

Im 

Luminous gain 

22,500 

42.000 

Equivalent 
background input 

1.8 X 10“'^ 
Im/cm^ 

1.04 X 10’“ 
lm/cm‘ 


could effectively be switched much faster if the voltages 
across the microchannel plate (MCP) and between the 
MCP and phosphor were left on, and only the voltage 
between the photocathode and the MCP was switched. 
This was accomplished with a 200-V negative pulse which 
was varied in width from 0.5 to 82 jusec. In other experi¬ 
ments at LASL, a similar tube has been successfully gated 
for as short a time as 5 nsec. This gating capability allows 
the radiographic sensor to operate in high background 
orbits, since a direct-current background (which appears 
as light on the screen) can easily be overcome if the image 
intensifier is operated for only several microseconds in 
synchronization with a beam burst from a flash x-ray 
machine or linac. 
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The present radiographic system does not use a 
mirror to reflect the image by 90°, and the image intensi* 
Tier is therefore not shielded from the beam. The system 
was used in this manner throughout the entire experimen¬ 
tal program without any apparent ill effects. In a flyable 
system the in-line camera has the distinct advantage of 
saving this extra shielding weight. 

All permanent radiographic data records were in the 
form of Polaroid pictures taken by direct film contact 
with the output fiber optics of the image intensifier tube. 
Much of the work was done with Polaroid 3000 film, but 
slower, high-contrast films were also used. Pictures for 
reproduction were taken on Polaroid Type 55 positive/ 
negative film. Switching from one film to another re¬ 
quired, of course, compensations in radiation dose, tube 
gain, etc. 

Comparisons of film sensitivity and vidicon sensitiv¬ 
ity were made to provide some basis for relating the 
performance of the system with that of a flight system 
with video recording. Several references cite a nominal 
vidicon sensitivity of 0.1 fc illumination on the vidicon 
target. A 30-Hz video scan corresponds to an exposure- 
per-frame of 3.3 x ICT^ fc^ec or 3.55 x ICT^ metercandle- 
sec. Using data sheets from the Polaroid Corporation, it is 
found that the Type 57 film (3000 ASA) requires 
3.98 X 10"^ metercandle-sec to produce a 25% contrast 
picture; and that the Type 55 positive/negative film 
(50 ASA) requires a 3.6 x 10“^ metercandle-sec for a 25% 
contrast picture. A 25% contrast picture is subjectively a 
good picture and was therefore chosen for the compari¬ 
son. The above figures show that the 3000-speed film is 
about 100 times more sensitive than the nominal vidicon, 
and that the posidve/negative film has about the same 
sensitivity as the nominal vidicon. An improvement in the 
vidicon gain by 100 is probably easily achievable in a 
flight system by slow scan or other techniques, so the 
sensitivity of the radiographic system with 3000-speed 
film is comparable to that of a flight system. 

II. Lens Evaluation 

Each of the different lenses to be evaluated present¬ 
ed a different image format and a different focal length. 
This latter dimension is critical in that the photocathode 
of the image intensifier is at high voltage and arcing can 
occur. The lens evaluation also provided observations of 
the behavior of a very wide-angle lens placed close to the 
screen. The five lenses in Table B-II were evaluated using 
^Co as the radiographic source. 

The f/0.78 lens was used with only panial success 
because the back focal length was small eno^ogh to cause 
intermittent arcing from the inuge intensifier. From the 
data that were taken, however, it was approximately 


TABLE B-II 

CHARACTERISTICS OF EVALUATED LENSES 



Focal 

View 

Lens-to- 



Length 

Angle 

Screen 


Aperture 

(mm) 

(deg) 

Disunce 

Make 

f/0.78 

25 

30 

6 ft 

Canon 

f/0.95 

25 

30 

6 ft 

Angenieux 

f/0.95 

50 

30 

6 ft 

Canon 

f/1.3 

15 

-^7 

45 in. 

Angenieux 

f/1.8 

5.9 

94 

18 in. 

Angenieux 


twice as fast as the f/0.95 lens. The 25-mm, f/0.95 lens 
was used nornully with the 18-mm image intensifier, 
because the image size from this lens closely matched the 
size of the intensifier. Similarly, the 50-mm, f/0.95 lens 
was used normally with the 23-mm image intensifier. 
However, during the lens comparisons, all the lenses were 
used with the same tube. The two f/0.95 lenses appeared 
to perform comparably except for the difference in image 
size. 

One concern with the wide-angle lenses is illustrated 
by the f/1.8, 94° case in which the lens was 18 in. from 
the center of the screen and 31 in. from the extreme 
edges. This means that the center of the screen should be 
contributing three times more light than the edges, and it 
was expected that the radiographs would have a bright 
area near the center. The effect, if present at all, was not 
apparent in any of the test radiographs made. In addition, 
any distortion due to the extreme close-up was not in 
evidence. These tests also indicated that the nominal f 
numbers of the lenses were apparently not affected by the 
severe geometries. 

III. Screen Evaluation 

Three screens were evaluated using a 600-mCi ^Co 
source and exposure times of 30 sec. The screens were 
0.010 in. and 0.020 in. ZnCdS (Radelin, PFG) and 
0.020 in. CaW 04 (Radelin STF-2), and several piaures 
were taken with each screen for equivalent exposure 
times. The 0.020-in. ZnCdS and 0.020-in. CaW 04 screens 
had twice the light output of the 0.010-in. ZnCdS screen, 
and it was found that the 0.020-in. ZnCdS screen had an 
extraordinarily long persistence when exposed to daylight 
or to radiation. The other two screens did not appear to 
suffer from this persistence problem. The thicker screen 
(0.020 in.) was made of two thinner layers bonded to¬ 
gether, a manufacturing technique which may have con¬ 
tributed to the high persistence. 
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CaW 04 screens have some advantages over the 
ZnCdS screens. They are more sensitive to harder radia¬ 
tion, and the blue light output (4400 A) is better mated 
to the S-20 and S-25 photocathodes of the image intensi- 
fiers, compared with the yellow-green (5500 A) of 
ZnCdS. Also the light output as a function of temperature 
is quite different for the two screens. These phenomena 
favor the CaW 04 screen for this application and will be 
discussed further. 

The optimum thickness of converter material 
(either lead, tungsten, or tantalum) in front of the scintil¬ 
lating screen is a strong function of the characteristics of 
the source used. For example, 0.040-in. lead sheets are 
used routinely at the Pantex linac, whereas 0.010-in. lead 
sheets are used routinely at LASL’s GMX-1 betatron. 
Thicker lead converters give fuzzier images, however 
thicker converters also give bener shielding to the scintil¬ 
lating screen (or film) in the presence of extraneous 
backgrounds (room scatter, target scaner, and source- 
associated background). 

Several thicknesses of lead, ranging from zero to 
0.060 in., were used with CaW 04 screens at the SCAPP at 
4 and 7 MeV. The quality of the radiographs improved 
rapidly up to 0.040 in. and slowly thereafter. However, 
the image obtained using 0.010 in. of lead was quite 
acceptable, whereas the radiograph taken with no lead 
Lacked contrast. Because the SCAPP was found to have a 
great deal of source-associated background, the thicker 
sheets of lead provided more shielding rather than acting 
as converters. 


IV. Screen-Lens-Photocathode Matching 

Optimization of an intensifier-typc radioactive 
system requires careful matching of the screen output 
spectrum with the spectral transmission characteristics of 
the objective lens and with the spectral response charac¬ 
teristics of the intensificr photocathode. Because most of 
the radiographic experiments were done with a non¬ 
optimum combination, the relative response of other 
combinations of these components was calculated to give 
some indication of the performance increase to be expect¬ 
ed. This information would certainly temper a judgment 
on such matters as optimum exposure levels or absolute 
system sensitivity for this type of radiographic system. 

The basic calculation of the response of a given 
intensifier photocathode and a given objective lens when 
exposed to light of some given spearum is 



where I is the photocurrent emitted by the photocathode, 
S(X) is the power spectrum of the incident light before its 
incidence on the lens, T(X) is the spectral transmission of 
the lens, and R(X) is the spectral response or absolute 
sensitivity of the photo cathode. The limits of integration 
are over the incident light spectrum. The absolute spectral 
radiant power being emined from the scintillating screen 
is not well known, but this does not prevent calculation 
of relative system outputs because the spectral shapes of 
CaW 04 and ZnCdS light emissions are well known. 

Figure B-1 shows the spectral shapes of CaW 04 and 
ZnCdS light emissions, as well as the variety of photo¬ 
cathode response functions selected for trial calculations. 
For all calculations a **typicar* lens transmission function 
was used and is also shown in Fig. B-1. This lens transmis¬ 
sion function is that published^ for a Canon 50-mm, 
f/0.95 lens, and is representative of fast, medium focal 
length lenses. The photocathode response functions 
selected included the S-20, the 116 which favored PFG, 
the 111 which favored CaW 04 , and the 135 which was 
flat across the specual band encompassing both CaW 04 
and ZnCdS. The number designations are RCA nomencla¬ 
ture. 



Fig. B-1. 

Spectral functions of screens, lens, and photo- 
cathodes. 
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A simple computer program was written to inte¬ 
grate the product of a spectral funaion, the transmission 
function, and a photocathode response function. This 
calculation was done for all combinations of the functions 
shown in Fig. B-1, and the results of these calculations are 
shown in Table B-III. Note that with the same lens there 
is a range of about a factor of two in the photo cathode 
output depending on the screen-photo cathode matching. 
This particular lens favors the ZnCdS by removing about 
20% of the CaW 04 emission at the spectral peak. 

For the present experiments note that for the same 
radiant power emission at the screen, CaW 04 outperforms 
ZnCdS by a factor of 0.887/0.642 or 1.38 in terms of 
measured response, i.e., photomultiplier response, for in¬ 
stance. The relative light production efficiency in the two 
screens has not been considered here. However, the 
CaW 04 . because of its high average atomic number is a 
more efficient light producer than ZnCdS for the energet¬ 
ic photons of interest in this program; so considerations 
of relative light-producing efficiencies will tend to en¬ 
hance further the performance capability of the CaW 04 
screen. The importance of this ratio is that a CaW 04 
screen with the same light output (watts per square centi¬ 
meter, for instance) as a PFG screen would provide 38% 
more response in our image intensifier which happens to 
have the S-20 photocathode. 

V. Cold Screen Experiments 

Experiments were performed to investigate the be¬ 
havior of CaW 04 and ZnCdS scintillating screens over a 
range of temperatures. Information from several sources 
indicates that the light output from these screens may 
vary drastically with temperature, a phenomenon of spe¬ 
cial interest for space applications of the materials. In 

TABLE B-III 

CALCULATED SCREEN-PHOTOCATHODE RESPONSE 


Relative Response 


Photocathode 

CaW04 

ZnCdS 

no (S20) 

0.887^^^ 

0.642^"^ 

111 

0.946 

1.000 

116 

0.992 

0.482 

135 

0.752 

0.663 


(Response)CaW04/S-20 

(R«Ponse)ZnCdS/S-20 


these experiments small samples of each screen material 
were cooled with liquid nitrogen while the screen was 
being irradiated with a ®®Co source and while the light 
output from the screen was being measured with a photo¬ 
multiplier and electrometer. The dose rate delivered to 
the screen for all trials was about 85 R/h. 

The design of the experimental geometry permitted 
(1) low thernul impedance between the screen material 
and the cryogen, (2) lighttightness to insure low optical 
backgrounds, and (3) insulation of the photomultiplier 
from the cold parts of the assembly. As the cryogen was 
introduced into a dewar containing the cold finger, both 
the output of the photomultiplier and the temperature of 
the cold plate were recorded. The initial temperature was 
either room temperature or as high as 50®C which was 
achieved by heating the screen-backing with a heat gun. 
When the temperature of the plate reached about -lOO^’c, 
heat was again applied to the plate with the heat gun. 
Light output vs screen temperature was again recorded 
during the warm-up of the assembly. Careful control 
measurements using artificial Light sources were made to 
insure that the photomultiplier response was not being 
affected by temperature excursions during the exper¬ 
imental trials. During such measurements the output of 
the photomultiplier was constant over the extremes of 
temperature. 

The results of these experiments are shown in 
Fig. B-2. The light output from CaW 04 increases with 
decreasing temperature, whereas the reverse is true for 
ZnCdS. The CaW 04 was well-behaved in that the light 
output tended to follow the same curve as the tempera¬ 
ture was lowered and raised. On the other hand, ZnCdS 
exhibited a peculiar hysteresis effect and the warming 
curve tended to jump to a new locus which was higher 
than the cooling curve. The behavior of the ZnCdS seem¬ 
ed to depend on recent thermal history; for instance, the 
initial light output at room temperature was maximum 
whenever the material had soaked at room temperature 
for many hours, but during subsequent temperature cy¬ 
cling this initial high output at room temperature was not 
achieved. It was also observed that the behavior of the 
material during each temperature cycle varied from the 
preceding cycle. Although the experiments were not con¬ 
tinued long enough to establish it unequivocally, it 
appeared that after several temperature cycles the mate¬ 
rial was approaching a well-defined characteristic curve 
which was single-valued for each temperature. As a con¬ 
sequence of this anomalous behavior, the characteristics 
of ZnCdS are displayed in Fig. B-2 as a broad band 
instead of a curve. 

To relate these two curves it should be noted that 
the experimental geometry, radiation source-screen spac¬ 
ing, was altered between trials to maintain the same signal 
output from the photomultiplier at room temperature. 
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Fig. B-2. 

Screen output vs temperature for CaWO^ and 
ZnCdS. 


about 25**C. This temperature is the approximate crossing 
point for the two curves, and provides the absolute com¬ 
parisons between the two materials. Recall, however, that 
in a real satellite inspection system the effects of the 
objective lens and the response of the specific photo¬ 
cathode used with the image intensifier must be taken 
into account; in these experiments the recording photo¬ 
multiplier had the S-20 photocathode. Further, it must be 
realized that different radiation absorption efficiencies at 
<iiffcrent photon energies will change the relative response 
of the two materials at any temperature. 

Light output variations with temperature are signifi¬ 
cant, about a factor of two over the temperature ex¬ 
tremes. The observed changes are probably not so gross as 
to eliminate either screen as a candidate material. How¬ 
ever, in the design of a flight system, these data point to 
thermal control of the radiographic screen as one way of 
enhancing system performance. 


VI. Summary 

The components of the radiographic camera, which 
i ncluded scintillating screen, converter, obje ctive lens, 
image intensifier, and power supply, have been thorough¬ 
ly demonstrated. Ei^er the 0.020-in. CaW 04 or the 
0.020-in. ZnCdS (if the persistence problem can be 
solved) screen is suitable, but the CaW 04 screen is some¬ 
what superior. The converter material should be strong 
enough to suppon the scintillating screen and still be 
flexible. Lead, being a soft material, has very little 
strength in thin sheets. Materials such as tungsten or 
tanulum have been used as converters and are much 
better suited to this application. Because tantalum and 
timgsten have higher density but lower atomic number 
than lead, an equal thickness can be used for all practical 
purposes. The optimum thickness of lead converter was 
found to be 0.040 in. (using the SCAPP) but radiographs 
of fair quality were made with 0.010-in. lead. Because the 
weight of the converter becomes appreciable for large 
screens, a 0.010-in. thickness of either tungsten or tanta¬ 
lum is recommended. 

Because the f/0.95 lens was used typically in the 
program, all of the recommended intensities, etc., are 
based on a lens with this speed. A faster (or slower) lens 
would require proportionately less (or more) radiation to 
give the same picture. The fastest lenses used (fA).95 and 
f/0.78) had also the narrowest view angle (30®), and the 
slowest lens (f/1.8) had the widest angle (94®). Assuming 
that these are the best lenses available, then the choice of 
lens becomes a trade-off between the intensity of the 
source and the distance from screen to lens. For example, 
if the screen size chosen is 6 by 6 ft, then the wide-angle 
lens is 3 ft from the screen and the narrow-angle lens is 12 
ft from the screen. Because it should be relatively easy to 
design a lightweight, but rigid, light baffle for a 12-ft 
screen-to-lens separation, the f/0.95 lens (or faster) is 
recommended. The f/0.78 lens (also 30® view angle) is 
twice as fast as the f/0.95 lens, but is not recommended 
because of the arcing problem. 

The proximity-focused (generation 2) image intensi- 
fier can easily be coupled directly to a vidicon with a fiber 
optics bundle rather than viewing the output with another 
lens, which costs another factor of at least 100 in inten¬ 
sity. It is also apparently unaffected by modest amounts 
of radiation, even at extremely high dose rates. As a result 
of this the lens and image intensifier can view the screen 
direaly and are therefore unshielded (except by the ur- 
get). 
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APPENDIX C 

SPACE RADIATION BACKGROUND 


I. Introduction 

Nuclear detectors arc sensitive to the primary and 
secondary effeas of a number of space radiation back¬ 
grounds. The most intense effects arise from energetic 
charged panicles trapped in the earth’s magnetic field. 
Occasionally intense solar-proton flares produce energetic 
protons of intensities approaching those of the trapped 
particles. Galactic cosmic rays, consisting mostly of high- 
energy protons, affect the detectors by deposition of 
charge in sensitive regions and by secondary-particle pro¬ 
duction in both the earth’s atmosphere and the space¬ 
craft. 

Neutrons which leak back into space from the at¬ 
mosphere are the only important source of neutrons in 
free space, while neutrons produced in the spacecraft 
nuterials add to this free space background. The last 
background to be considered, that of cosmic gamma rays, 
is low enough to be almost insignificant. In the following 
sections the backgrounds will be treated in increasing 
order of importance. 


II. Cosmic Gamma Rays 

Cosmic gamma rays above 0.25 MeV have been 
measured by Vene et al,^’ who measured a flux of 0.7 
photons cm'*sec"' above 0^5 MeV and 0.2 photons 
cm'* sec'' above 0.6 MeV. Above 0.6 MeV the intensity 
falls as . Because the flux is constant, except for 
variation of eanh shadowing with satellite altitude, sub¬ 
traction of this background should be simple if it becomes 
necessary. 


III. Solar Proton Flares 

Solar particles of energies from a few million elec¬ 
tron volts up to a few gigaelectron volts arc associated 
with certain solar flares. These particles arc predominant¬ 
ly protons with varying fractions of alpha particles and 
heavier nuclei and arc also called solar cosmic rays. The 
large events, occurring a few times a year, can produce 
fluxes exceeding the galactic cosmic rays by many thou¬ 
sands of times and can last one or two days. During such 
large events, neutron- and*gamma-ray detectors would be 
unusable except in the low magnetic latitude portion of 
the lower orbits. There the earth’s magnetic field would 


shield out all except the most energetic particles, with the 
penetrating flux depending strongly on the energy spec¬ 
trum of the particles. Thus, solar protons would negate 
neutron- and gamma-ray detection only rarely and prob¬ 
ably not at all for the midlatitudcs of the 300-nm and 
lower orbits. 


IV. Galactic Cosmic Rays 

At the top of the atmosphere, protons contribute 
almost 90% of the cosmic-ray flux. Helium nuclei contri¬ 
bute 9%, with heavier nuclei and electrons each contribut¬ 
ing 1%. The mean energy of the spectrum is about 1 GeV. 
The interstellar flux is constant and isotropic, but at the 
earth the flux is affected by solar aaivity. Short period 
decreases (Forbush decreases) of up to 15% in neutron 
production by the primary particles can occur after geo¬ 
magnetic storms which arc caused by solar activity. The 
flux gradually returns to the initial value in a day or two. 
Longer period variations arc associated with the 11-yr 
solar cycle. The flux decreases about 40% from the solar 
activity minimum to the maximum. Probably the increase 
in particles leaving the sun at solar maximum increases the 
trapped magnetic field carried along by particles. Thus, 
the magnetic field in interplanetary space increases and 
excludes more of the cosmic rays, thereby causing the 
decrease in flux. 

Consider first the atmospheric interaction. The pri¬ 
mary particles produce protons and neutrons and other 
particles by nuclear reactions with the air molecules. The 
protons lose energy more rapidly to the air because they 
are charged. The neutrons are scattered, with a resultant 
energy loss, and some neutrons leak back out the top of 
the atmosphere. The resultant neutrons above the atmo¬ 
sphere are called albedo neutrons. The number of primary 
cosmic-ray particles reaching the top of the atmosphere at 
any point is affected by the earth’s magnetic field. A 
relatively low-energy particle can reach the earth at the 
magnetic poles, but particles must *have much higher ener¬ 
gy to penetrate to the earth at the magnetic equator. 
Consequently, the cosmic-ray flux at the top of the atmo¬ 
sphere is at a minimum at the magnetic equator and rises 
gradually to a maximum in the magnetic polar regions. 
The albedo neutron flux produced by these particles 
varies by a factor of 7.4 between the equator and the 
poles.^ The neutron flux at the top of the atmosphere 
(50 km) was determined to be 0.52 cm“^sec*‘ at the poles 
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and 0.070 cnT^scc"' at the equator at solar maximum. 
The flux variation with altitude was consistent with an 
RT^ dependence, where R is the distance from the center 
of the earth. This dependence is what one would calculate 
for an isotropic leakage flux. 

The effeas of cosmic rays on the detector- 
spacecraft system are difficult to calculate, because the 
geometry is complex and the high-energy particle- 
production cross sections are not well known. The effects 
arc much more significant for the neutron detector than 
for the other detectors, so some experimental results for 
the Vela satellite neutron detectors will be used. These 
detectors arc moderated ^ He proportional counters with 
no guard counters. The moderators arc 6.5-in.-diam poly¬ 
ethylene cylinders, 6.5 in. long. The efficiency for fission- 
spectrum neutrons is 10 cm^ (about 4%). The Vela orbits 
are circular at approximately 18 earth radii, and albedo 
neutrons arc significant at this altitude. The cosmic-ray 
flux is approxinutcly that of the lower altitude polar 
regions since the earth’s magnetic field shielding is insig¬ 
nificant. The detector count rate from all cosmic-ray- 
produced effects varies from 8 sec”* at solar maximum to 
14 sec”* at solar minimum; equivalent to a neutron flux 
of 0.8 to 1.4 cmT^scc”*. The spacecraft weight is approxi¬ 
mately 350 kg, and the detector is well inside the space¬ 
craft. At the top of the atmosphere, this detector- 
spacecraft combination would produce a count rate equiv¬ 
alent to a neutron flux of 1.3 cm”^scc’‘ (0.8 +0.52) at 
the poles and 0.17 cm”^scc”‘, 

0.8 + 0.52 \ 

7.4 ) * 

at the equator at solar maximum and 5/3 of these rates at 
solar minimum. The local effects produced in the detec¬ 
tor-spacecraft system are assumed to vary directly with 
the albedo neutron flux. This follows from the assump¬ 
tion that cosmic-ray neutron production in spacecraft 
materials is proportional to neutron production in the 
atmosphere. 

V. Trapped Particles 

The Van Allen or radiation belts of the earth are 
regions containing intense fluxes of energetic protons or 
electrons, or both. These particles are trapped by the 
earth’s magnetic field, which is approximately a dipole 
field. Particles execute spirals along magnetic field lines 
which converge in the magnetic polar regions. The angle 
between the velocity vector and the field line (the pitch 
angle) increases as the particle continues in the direction 
of increasing field until it reaches 90®. At this point (the 
mirror point), the spiral motion along the field line 


changes direction by 180® and continues until the particle 
mirrors in the other hemisphere. If a mixture of particles 
is inserted on a field line at the magnetic equator (the 
point of lowest magnetic field intensity for any field line), 
the particles will spiral along this line and mirror at 
various points. Those with smaller pitch angles will mirror 
closer to the poles than those with larger pitch angles. The 
mirror points are independent of energy. Thus, the parti¬ 
cle flux on this line would be greatest at the equator and 
would decrease as the field increases toward the poles. 
For any field line the particle flux should be a function of 
the magnetic field only. 

Consideration of all identical field lines (identical 
except for their longitude) gives a shell which interseas 
the eanh in a circle about each magnaic pole and is 
farthest from the earth at the magnaic equator. The 
distance in earth radii from the center of the earth to the 
shell equator is called the L value of the shell. Particles 
trapped on any field line of a particular L shell drift 
longitudinally around the eanh in paiods of the order of 
1000 sec while mirroring between the hemisphaes in 
about 1 sec. The trapped particles for a particular L shell 
mix thoroughly with all others in that shell. Thus, the 
trapped particle flux can be charaaerized by two vari¬ 
ables: the magnaic field, B, and shell value, L. 

Vette and others^* have correlated data from many 
satellites which penetrate the radiation belts and have 
devised models of particle environment in terms of B and 
L. For each L one additional piece of information is 
needed, namely, the particle energy speara. These are 
obtained experimentally for each shell. The NASA Na¬ 
tional Space Data Center has incorporated model environ¬ 
ments on punched cards and has written a program to 
calculate satellite orbits, to daamine B and L from 
position, and then to daamine electron and proton 
fluxes and speara from the model environments. The 
program is called TRECO. 

The TRECO code was modified to record spectral 
information about the model radiation environment at 
seleaed time intervals along an orbit. In the case of 
trapped protons, the spectral information was obtained 
by plotting on film the particle flux in selected energy 
intervals. For the elearon belts, funha analysis was re¬ 
quired so the complete spearum at seleaed places was 
also punched on cards for later use. 

The model trapped radiation environments which 
were used were also obtained from NASA and are referred 
to as the Vene models AP-1, AP-6, AP-7, and AE-2. The 
models are static and therefore, in the particular case of 
the electron model AE-2 of 1964, the model does not at 
the present time reflea the avaage composition of the 
belts. Variations in the avaage composition of the elec¬ 
tron belts, which are expeaed to have taken place since 
the model year 1964, are a decrease in inna belt fluxes 
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caused by decay of the Starfish electrons and an increase 
in outer belt fluxes accompanying the increase of solar 
activity in the solar cycle. The inner belt discrepancy is 
the only one about which much information is available, 
so the AE-2 nK>dei environment was modified to approxi- 
nute the present situation by truncating each spectrum 
for L < 2.5 above the energy 0.75 MeV. This corresponds 
to the NASA model electron environment which is 
expected to be published in 1972. 

VI. Gamma Detector Background 

Further analysis of the electron data was necessary 
because detector background is directly produced through 
bremsstrahlung photons. To calculate the bremsstrahlung 
radiation, an empirical method^^ was used in which the 
measured bremsstrahlung intensity produced in a thick 
target in the forward hemisphere is related to the ratio of 
the energy of the photon produced to the energy of the 
incident electron. Using a single empirical curve, a trans* 
mined bremsstrahlung spectrum can be calculated from 
the incident electron spectrum. The method was checked 
using the Grumman electron spectrum^ and did give the 
same photon spectrum that the Grumman report con- 
uins. To calculate a bremsstrahlung spectrum in the back 
hemisphere, a backscaner curve was obtained by taking 
the difference between published 47r and forward hemi¬ 
sphere curves. 

To calculate the amount of bremsstrahlung radia¬ 
tion transmined through shielding nuterials in a simple 
manner, some provision must be nude for radiation which 
is initially scattered but is eventually counted because of 
the spatial extent of the active detector volume or of 
multiple scattering. Such scattering was accounted for by 
using as an absorption cross section the total energy 
absorption cross section plus one quaner of the total 
scattering cross section. This procedure has been shown to 
be valid for radiation in the range 10 to 100 keV and, for 
lack of other means, has been applied between 100 keV 
and 7 MeV. 

The specific geometry calculated corresponds to a j 
shielded Gc(Li) detector. The outer layer of 2-in. LiH 1 
serves to stop the incident electrons while producing a ] 
minimum amount of bremsstrahlung because of the low 
average Z. Next, a 1/4-in. layer of lead attenuates the 
large number of low-energy photons. The guard counter 
follows and is 2-in.-thick Csl. The remainder consists of a j 
3-in.-diam, thin Mylar jacket enclosing the 45-mm-diam 
by 50-mm-long active volume. The length of the shield 
was 20 in. and the field of view of the collimator 10®. The j 
collimator aperture was closed to direct entry of charged 
particles by 2 in. of LiH which was followed by a soft 


photon filter composed of 1/2 g/cm^ of molybdenum and 
1/6 g/cm^ of tin. 

The assumption was made that the entire detector 
aperture would view the target, making the background 
calculation independent of distance. It was also assumed 
that charged particles would not enter the collimator 
during inspection. This could be accomplished in the 
lower orbits by not looking in directions approximately 
perpendicular to the magnetic field since the charged 
particle pitch angles are close to 90®. Thus, the back¬ 
ground entering through the collimator was assumed to be 
back hemisphere bremsstrahlung produced in the near 
side of a thick aluminum target. The cylindrical sides of 
the shield were assumed to be directly exposed to the 
charged-particle environment, and the ends of the shield 
were assumed to be far enough away to produce negligible 
background. 

During the calculation, all of the energy absorbed in 
the guard counter is tallied to produce the guard counting 
rate, and the Gc(Li) background rate is found from the 
photon spectrum reaching the active volume using a 
sample 6% Ge(Li) intrinsic photopeak efficiency curve. 
One must realize that the calculations of photon number 
and energy incident on the Ge(Li) have been simplified to 
permit calculations of a large number of cases. An accu¬ 
rate calculation would require an elaborate photon trans¬ 
port code (Monte Carlo) configured to a particular de¬ 
tector. Therefore one should use the heavily-shielded de¬ 
tector count rates as order-of-magnitude estimates. 
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APPENDIX D 

NEUTRON SPECTRUM UNFOLDING 


The use of gamma-ray interrogation either actively 
or passively allows the identification of elements with low 
(7,n) reaction thresholds by the deteaion of the resultant 
neutrons. The elements of interest in the present applica¬ 
tion are and ^Be which have (7,n) thresholds at 2.226 
and 1.666 MeV, respectively. The elements can be distin¬ 
guished by the maximum energy (endpoint) neutrons 
produced. These energies are given below. 

2.226) 

“ 8/9(E.y - 1.666) 

For a simple system a well calibrated (in energy/channel) 
pulse-height analyzer would be adequate to distinguish 
these endpoints. 

More information could possibly be obtained by 
unfolding the pulse-height spectra from the detector to 
obtain the neutron spectra incident on the detector. It 
may simplify identification in the case of significant 
shielding or moderation and in the presence of high- 
energy fission neutrons. In addition, it might be possible 
to obtain more quantitative infornution from the number 
of neutrons that can be assigned to each component 
[^Be( 7 ,n), ^D(7,n) and fission]. To this end the experi¬ 
mental and numerical techniques involved in unfolding 
pulse-height spectra obtained from ^He proportional 
counters were investigated and are described here. 


The equation for the response of a detector, D(?), 
to a signal, 0(E), can be written as 

D({) = f A(E,J)0(E)dE , (D-1) 

0 

where A(E,f) is the response at ^ to a unit signal at E. In 
the present application (neutron-spectrum measurement) 
0(E) is the neutron spectrum as a function of energy E, 
and D(|) is the pulse-height spectrum as a function of 
pulse height f. This equation is commonly solved by 
dividing both the range of E and the range of $ into a 
finite number of neutron-energy and pulse-height incre¬ 
ments to obtain a set of simultaneous linear equations: 

J 

Di= L Aij«j (D-2) 

or 

A0 = D 

Anempts at solving Eq. (D-2) by inverting the matrix A 
( 05 = A”* D) may lead to unbounded errors in the solu¬ 
tion 0^. These are a result of the nature of the equation 
and the fact that D is not known exactly but only as 

Di+aj. 

The least-structure method solves Eq. (D-2) subject 
to the constraints (l)that the (2nd derivative)^ of the 
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solution be minimized and (2) that the residual e [defined 
by Eq. (D-3)l be consistent with the 

e*A0j-D (0-3) 

error with which D is measured. Consistency here will 
mean 



The second derivative will be expressed by the following 
numerical approximation to the second derivative: 





The constrained minimization is accomplished by intro¬ 
ducing the Lagrange multiplier X and minimizing the 
function: 



or solving the set of equations 


30 

30) 


= 0 


forj = 1 ... J 


It is assumed that 


3 ^ 0 3^0 


Forcing the solution to have a minimum second 
derivative as well as to satisfy Eq. (D-2) tends to prevent 
the unbounded oscillations mentioned above. The amount 
of smoothing used in the solution (how small the second 
derivative is forced) is strongly dependent on the estimate 
of the errors, (7, in the measured data, D. Clearly, too 
much smoothing will eliminate actual structure in the 
solution. 

As a test of these techniques for obtaining neutron 
spectral information from pulse-height distribution meas* 



NEUTTON ENERGY (MaV) 

Fig. D-l. 

Fission neutron spectrum from ^Pu multiplying 
sphere. 

urements using ^He proportional detectors, the neutron 
spectra from a ^’Pu sphere, a ^^Cf source, and from the 
(7,n) reaction on ^Be and at two electron energies 
were obtained. These measurements, shown in Figs. D-1 
through D-3 and Fig. 43 were intended only as a proof of 
the experimental and numerical techniques and not as a 
demonstration of the full capability of the technique. It is 
felt that more can be done with the use of rise-time 
discrimination, with better measurements of the detector 
response function, A(E,{), and by using response func¬ 
tions more closely spaced in energy to take advantage of 
the resolution of the detectors. 

In addition, by obtaining normalized response func¬ 
tions (counts/channel/incident neutron), absolute spec¬ 
trum measurements can be made. The spectra presented 
here have arbitrary normalizations. 
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NEUTRON ENEROY (M«V| 

Fig. D-3, 

Deuterium and beryllium pbotoneutron spectra, 
4.5-MeV bremsstrablung. 


* 0.1 1.0 to 

NEUTRON EN^GY (MtV) 

Fig. D-2. 

Spontaneous fission neutron spectrum, ^^Cf. 




APPENDIX E 

MONTE CARLO ANALYSIS OF *He NEUTRON DETECTORS 


I. Introduction 

General response characteristics of ^ He neutron de¬ 
tectors were studied with Monte Carlo techniques. These 
studies complement the experimental evaluation and ap¬ 
plication of ^He detectors for neutron counting and spec¬ 
troscopy. Effects of detector geometry, fill-gas composi¬ 
tion, and gas pressure on detector response were exam¬ 
ined by parametric Monte Carlo calculations. The calcu¬ 
lated detector responses provide a guide to optimizing 
detector design for the particular applications considered 


in this study. A discussion of calculational techniques and 
typical results follows. 

II. Generation of Range-Energy Table 

A range-energy table is a basic requirement for the 
Monte Carlo calculation of recoil spectra. For a ^He 
detector, this table of (E,R) values gives the range R(E) of 
the recoil alpha particle in ^He, where E is the initial 
energy of the alpha particle. The range is calculated from 
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where dE/dx is the alpha-particle stopping power of ^He 
at the appropriate gas pressure. Values of stopping power 
were obtained from the recent and very useful compila¬ 
tion of Northcliffe and Schilling.^^ With the format of 
Ref. 33 as a guide, R(E) was evaluated analytically in the 
range O < E < 0.05 MeV using -dE/dx = kE'^, with k de¬ 
termined from the stopping power at 0.05 MeV. For the 
energy range 0.05 < E < 48 MeV, R(E) was numerically 
evaluated from Eq. (E-1) using a seven-point Newton- 
Cotes formula to integrate tabulated (dE/dx)"* data. To 
provide a finer logarithmic energy scale, R(E) was evalu¬ 
ated at the intermediate values 



Fig. E-L 

Monte Carlo calculational geometry for isotropic 
point source of neutrons located on the axis of a 
cylindrical proportional counter. 

collision C to the nearest detector wall is com¬ 
puted. If r > R(E^), one event is tallied for the 
energy group AE containing If r < R(E^), the 
distance r is expressed as 


E = exp [0.5 (In Ej + In Ej^j)l , (E-2) 

where Ej and Ej ^ ^ are consecutive values on the refer¬ 
enced energy scale. 

111. Monte Carlo Procedure 

Monte Carlo calculations were made for an isotropic 
point source of neutrons located on the axis of a cylindri¬ 
cal proportional counter. The calculational geometry is 
shown in Fig. E-1. The Monte Carlo calculational steps are 
outlined: 

(1) Neutrons are emitted isotropically at S. Only 
those neutrons emitted within a source cone 
defined by the half-angle S are sampled. These 
neutrons enter the detector where they can es¬ 
cape from the detector with no collision, be 
captured in the detector wall, or be elastically 
scattered by ^He. A neutron having an elastic 
scattering collision at C will have an energy and 
direction (E„,Uj,,Vn,w„) after collision, as deter¬ 
mined from the differential scattering cross sec¬ 
tion d^o/dEdl2 (u„,Vj^, and w^j arc the direction 
cosines). 

(2) The energy and direction of the recoil nucleus 
(E^,Ujj,v^,w^) arc calculated from conservation 
of energy and momentum after collision. 

(3) Using a binary search interpolation routine, the 
recoil nucleus range R(E^) is computed from the 
rangc-cncigy table. 

(4) Step 4 treats wall and end effects. The distance,r, 
in the direction (u^,v^,w^), from the point of 


r* R(Ec^)-R(E^) . (E-3) 

where = energy of the recoil nucleus on 
reaching the detector wall and R(E^) » residual 
range, i.e., range of the recoil nucleus having 
energy E^. Since r and R(E^ arc known, the 
residual range is given by R(E^) » R(E^ - r. 
Using inverse interpolation, E^ is computed 
from the range-energy table and one event is 
tallied for the energy group E containing E^-E^. 

(5) Steps 1 through 4 arc repeated until enough 
neutron histories have been computed to give the 
desired statistical accuracy. The resulting tally 
per energy group gives a histogram of the ideal 
(perfect resolution) recoil pulse-height spectrum. 

Actiul recoil pulses produced in the proportional 
counter by monoencrgetic neutrons will have an energy 
distribution due to the statistical nature of the ion pro¬ 
duction and energy loss processes. Effcas of statistical 
processes, local field inhomogeneities, pulse rise-time, and 
external electronics on detector resolution arc difficult to 
calculate accurately from first principles. Rather than 
attempting to include these effects analytically, we treat¬ 
ed detector resolution by distributing the pulse height in 
each energy group with a Gaussian function. Detector 
resolution expressed as a Gaussian distribution has been 
shown to be a valid technique in the interpretation of 
measured ^Hc recoil pulse-height spectra. 

After completion of Step 5 above, the recoil pulse- 
height distribution is given by Cj counts accumulated in 
energy group AEj, j=l, —, J. Each pulse height is to be 
spread in energy by the Gaussian distribution 

Pj(E)*(l/v/2ira)cxp {-(E-Ej)^/2a^l . (E-+) 
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where Ej is the mean energy of group AEj and a is the 
standard deviation. The common measure of resolution, 
full width at half maximum, is E^ » 2.35a. Using this 
relation, Eq. (E-4) becomes 

Pj(E) » (0.938 /Ev 4) exp (-2.76 (E-Ej)VEv 4 ^J. (E-5) 

The distributed pulse height in group AEj, resulting from 
the computed ideal pulse height in group AEj, is calcu* 
lated numerically using Eq. (E-5) as 

Ei 

qj «y q Pj (E) dE, AEj = Ej - Ej.i . (E-6) 

Ei-l 

The total pulse height in group AEj is then given by 

J 

q* L Cij . (E-7) 

j-1 


amplitudes are ploned in Monte Carlo units of weight per 
million electron volt which, as noted, are proportional to 
counts/channel. In ^Hc, the maximum neutron energy 
transfer corresponds to 180** scattering, resulting in a 
recoil ^He nucleus having 64% of the incident neutron 
energy. 

As an example of the type of survey calculations 
that were made with these Monte Carlo techniques, 
effects of detector length on recoil pulse shape and ampli* 
tude were calculated. The response to 1.0-MeV neutrons 
for detector lengths from 5 to 25 in. is shown in Fig. E-3. 
In each case the source is located on the detector axis at 
1 m from the nearest detector end face. Figure E-3 shows 
that an increase in the detector length increases the peak 
amplitude at slightly less than a linear rate. The increased 
amplitude is due primarily to increased path length for 
maximum energy transfer (180® scanering). Other effects 
of detector length are summarized in Table E-I. 

The Table E-I data show that a reduction in de¬ 
tector length should not degrade the peak shape. Because 


Values of C^ from Eq. (E-7) give the recoil spectrum with 
detector resolution effects simulated by a Gaussian distri¬ 
bution. The resolution used in Eq. (E-5) is given by 

« 0.1 Ej. Ej < 1.0 MeV 

* 0.10 MeV, Ej > 1.0 MeV (E-8) 

This empirical relation is based on examination of typical 
resolution values reported in the literature. 

Figure E-2 shows alpha-particle recoil spectra result¬ 
ing from an isotropic point source of neutrons located 
1.0 m from the nearest detector end face. The deteaor 
length is 24 in., the o.d. is 2 in., and the ^He gas pressure 
is 8 atm. 

These calculations are based on the same neutron 
source strength, so that all recoil spectra in Fig. E-2 can 
be compared on the same relative amplitude scale. The 



Calculated recoil pulse-height spectrum in 8-atm 
^He detector for 0.4- to Z.O-MeV neutrons. 



Fig. E-3, 

Effect of detector length on recoil pulse spectrum 
for l.O-MeV neutrons. ^He detector, 2-in. o.d., 
8-atm pressure. 
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TABLE E-I 

EFFECTS OF DETECTOR LENGTH 


Detector 

Length 

(in.) 

Peak/V alley 

FWHM 

(MeV) 

Peak Location 
(MeV) 

5 

3.0 

0.20 

0.59 

10 

3.0 

0.20 

0.60 

15 

3.0 

0.19 

0.60 

20 

3.0 

0.19 

0.60 

25 

3.1 

020 

0.60 

only the 

number of counts 

in the peak 

is significantly 


changed by reducing the tube length, a short ^He detector 
could be used whenever the counting time is sufficient to 
achieve the required statistical confidence. 

The ^He detectors recommended for bremsstrah* 
lung-induced (7.n) and (7,0 interrogation are to have 
16-atm gas pressure. Monte Carlo calculated response 
functions for 1 b-atm detectors are shown in Figs. E-4 and 
E-5, for 1- and 3-MeV neutron sources, with tube of 2-in. 
o.d., length of 24 in., and source-to-detector distance of 
1 m. Increase in the gas pressure from 8 to 16 atm in¬ 
creases the amplitude of the full energy peak by a faaor 



Fig. E-4. 

Calculated receil pulse-beigbt spectra in 8- and 16- 
atm ^He detectors for l.OMeV neutrons. 



Fig. E-5. 

Calculated recoil pulse-height spectra in 8- and 16- 
atm ^He detectors for 3.0-MeV neutrons. 


of 1.5. The FWHM and peak/valley ratio are unchanged 
by a doubling of the ^He pressure. 

Calculated response functions discussed thus far are 
based on a source-to-detector distance of 1 m. To deter¬ 
mine what changes will occur when the detectors are 
moved to the baseline distance of 20 ft from the source, 
response functions were calculated for this source-to- 
detector distance. At the larger distance, the solid angle 
subtended by the detector end face is reduced by 0.027. 
The ratio of calculated peak amplitudes is 0.033. Thus, 
the intrinsic efficiency of the ^He detector is increased by 
0.033/0.027 = 1.22 in moving from 1 m to 20 ft away 
from a point source. This is because neutrons from the 
larger distance are entering the detector more nearly par¬ 
allel to the axis. The Monte Carlo calculations indicate 
that there should be no inherent problems in using 16-atm 
^He detectors at interrogation distances in the 
neighborhood of 20 ft. 

Detector responses were calculated for several neu¬ 
tron spectra and compared with measured pulse-height 
spectra from the same sources. Figures E-6 and E-7 show 
calculated ^He detector response to 4.5-MeV bremsstrah- 
lung-induced (7,n) spectra from D(D 2 0) and beryllium. 
The 8-atm, 2-in.-o.d.,^24-in.-long detectors were used with 
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Fig. E-6. 

^He detector response to 4.5-MeV electron 
bremsstrablung’induced (y^n) spectrum from D to 
D2O. 

a 1-m sourcc-to-dctcctor distance in the calculations and 
measurements. The computed detector response for 



Fig. E-7. 

^He detector response to 4.5-MeV electron 
bremsstrablung-induced (y,n) spectrum from beryl¬ 
lium. 

deuterium gives a good fit to the measured data. Some 
disagreement between the measured and calculated re¬ 
sponse is evident for the beryllium data. This difference is 
partially anributable to uncertainties in the contribution 
by excited state beryllium transitions to the photonuclear 
spectrum used in the Monte Carlo calculation. The general 
agreement between measured and calculated spectral re¬ 
sponse supports the application of photonuclear interro¬ 
gation to distinguish between beryllium and deuterium. 


APPENDIX F 

RADIOGRAPHIC IMAGE ENHANCEMENT 


Existing LASL facilities were used in an attempt to 
demonstrate the potential of image processing in restoring 
and extracting information from poor quality radiographs 
such as those telemetered to earth after a satellite inspec¬ 
tion. Image processing of television pictures taken by the 
Ranger and Mariner spacecrafts has already demon¬ 
strated the utility of the approach in a space application. 


Briefly, image processing techniques can be used to: 

(1) Correct for photometric distortion of the vidi- 
con, 

(2) Remove scan noise associated with slight varia¬ 
tions in the scanning beam position from scan to 




scan. 
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(3) Remove scan noise assodaced with the Gaussian 
shape of the intensity across the scanning beam 
dianKter, 

(4) Filter out recognizable, periodic system noise 
from other electronic subsystems on board the 
spacecraft, 

(5) Stretch contrast, 

(6) Sharpen edges for fine detail enhancement by a 
variety of schemes, 

(7) Smooth out white or random noise, and 

(8) Provide optimal (Wiener) filtering if adequate 
models of the noise and of the image proper can 
be made available. 

The above list is by no means exhaustive. Threre are many 
other processing techniques which might be applicable 
and useful. 

In this program photographic film instead of a video 
system was the basic image-recording medium; and since a 
TV camera is probably the primary noise source in a real 
system, image processing of film may not be especially 
applicable to the real system case. Nevertheless, it was felt 
that the potential of image-processing techniques to the 
enhancement of radiographs of real weapons might be 
demonstrable and that the opportunity to operate on 
such images should not be passed. 

Detailed discussions of image processing are to be 
found in Refs. 36 through 39, and the following discus¬ 
sion is meant merely to acquaint the reader with rudimen¬ 
tary notions of the image-processing techniques applied in 
this effort. 

The mechanics of image processing are straightfor¬ 
ward. An image, usually on a negative, is scanned by a 
rapid-scan microdensitometer which digitizes the image 
and stores it on magnetic tape. Computer storage limita¬ 
tions restricted the processing to 500 scans, each scan 
containing 500 discrete values. The stored image on tape 
then corresponded to a 500 x 500 array of density values, 
or 2.5 X 10^ density points. The 50^ aperture used for 
the scanning took full advantage of the inherent resolu¬ 
tion contained in the initial inuges. This inherent resolu¬ 
tion is determined by the image intensifier resolution and 
not the film, which has a much higher intrinsic resolution. 
The digitized image, now on tape, can then be computer- 
manipulated and the processed image “played back’* on 
Polaroid film. At LASL the scan and playback are per¬ 
formed on the same machine, a Photometric Data 
Systems Micro densitometer system 1050A, which con¬ 
tains a PDP-8 computer. 


Computer processing is typically done in either of 
two ways: (l)the digitized image (point locations and 
corresponding density values) is manipulated directly by 
the program, or (2) the density points are treated as a 
funaion which is represented by the points, and the 
function is transformed into the Fourier or spatial fre¬ 
quency domain before manipulation. 

With the first approach, for example, noise smooth¬ 
ing is possible by averaging each density point with several 
of its neighboring points. Contrast stretching is also done 
by working directly with the density values; the mean 
density value for the entire image is computed and then 
all density values on the tape are increased or decreased 
with respect to the mean value. There are many other 
processing schemes which involve the manipulation of 
only density data. 

The density data are transformed into a spatial 
frequency representation when some kind of filtering is 
preferred. Since the “function** to be transformed is rep¬ 
resented by discrete density values, the discrete Fourier 
transform (DFT) is used. One form of the DFT is 

N-1 

^ X|^ exp(- 27rjrk/N) r * 0, . . ., N - 1 , 

k*0 

where A^ is the rth coefficient of the DFT, X|^ denotes 
the kth density sample in a one-dimensional series, and 
j =\/nr. The A^*s are almost always complex, and the 
index r is sometimes called the frequency of the DFT. 

Operating on a series (a scan) of density values using 
the above expression produces a set of complex coeffi¬ 
cients representing the amplitudes of all the spatial fre¬ 
quency components which make up the “image” formed 
by the density values. From these calculated amplitudes, 
the power spectrum of an image can be generated. The 
power spectrum is nothing more than the square of the 
absolute value of the Fourier transform, i.e., 10^(0 P, 
where Dj(f) is the transform of d 5 (x), the density func¬ 
tion. The arguments f and x refer to frequency and 
[>osition respectively, and the units are conventionally in 
cycles/millimeter and millimeters. The power spectrum 
indicates the frequency content of an image and is useful 
in choosing the filter to be applied to the image in the 
spatial frequency domain. To relate the density values in 
real space and the power spectrum in the spatial fre¬ 
quency domain, it might be pointed out that the spacing 
of the density values turns out to be the wavelength of 
the highest frequency component in the resulting power 
spectrum; while the span of the density data is the wave¬ 
length of the zero frequency or fundamental mode in the 
power spectrum. 

The final detail of an image is contained in the 
high-frequency components of its power spectrum (once 
the image has been transformed). Similarly the broader 
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features and infomution about the overall exposure and 
contrast of an image are contained in the lower-frequency 
components of its power spectrum. This observation leads 
naturally to an enhancement scheme in which the image 
in spatial frequency domain is filtered in such a way that 
the high frequencies of the power spectrum are boosted 
and the low frequencies attenuated. The result, when the 
image is transformed into its original form, should be an 
enhancement of edges and of fine detail. 

Similarly certain kinds of periodic, predictable noise 
patterns may be evident in a given image, and may be 
easily filtered out by applying an attenuating filter in 
frequency space to that band of frequencies which 
represents the noise. If the filter has been properly chosen 
and properly applied, then transforming the image to its 
original form produces an enhanced inuge in which the 
bothersome noise has been removed. Such an approach is 
generally valid only when the power spectrum of the 
noise is in well-defined bands and does not coincide with 
the frequency bands in which essential information is 
contained. Otherwise the desired information will be fil¬ 
tered out along with the unwanted patterns. 

In this program, it was decided to use high pass- 
filtering and contrast stretching as the basic processing 
schemes. While it was not felt that these two schemes 
necessarily represent the optimum processing schemes, 
computer programs were readily available for them. Also, 
these techniques are well established and known to be 
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It is anticipated that effort on image enhancement 
will continue at LASL. Many different enhancement 
schemes will be attempted on the basic set of dau already 
coUeaed for this program. In addition, it is anticipated 


tl' c a video system will be assembled and video tape 
output, which more closely matches the data from a 
spaceborne inspector, will be available for the simulation 
ol noise suppression and information recovery. 
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